Glassy Slowing of Stripe Modulation in (La,Eu,Nd)2 a;(Sr,Ba)a;Cu04: 
A 63Cu and ^^^La NQR Study Down to 350 mK 
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®^Cu and ^"^^La nuclear quadrupole resonance and Zeeman perturbed nuclear magnetic reso- 
nance experiments are performed on the striped phase of the high temperature superconductors 
'La.2-x^s.xC\iOi and La2-a;-j/(Nd,Eu)i,Sra;Cu04. The first goal of the present study is to utilize the 
fact that ordered Cu magnetic moments exert a static hyperfine field on the ®^Cu and ^"^^La nucleii 
to deduce the charge density and ordered moment within the Cu02 planes. A hyperfine broad- 
ened NQR lineshape is observed in both La2-a:Bai:Cu04 and Lai.go-2!Euo.2oSra;Cu04 for a; ~ |. 
Detailed numerical analysis of the '''^Cu NQR lineshape establishes that widely accepted models of 
periodic sinusoidal or square- well shaped modulations of spin density waves with maximum moment 
~ 0.3^B, as inferred from elastic neutron scattering and /^SR measurements, can not account for the 
NQR lineshape unless we assume a relatively small ordered moment ~0.15/^s with a comparably 
large distribution. The second goal of the present work is to establish the temperature dependence 
of the fluctuation frequency scale of stripes. We find that the the fraction of missing ®'^Cu NQR 
intensity below charge ordering temperature Tcharge accurately tracks the temperature dependence 
of the charge order parameter as measured by scattering methods. By fitting a single model to the 
temperature dependences of the wipeout fraction F{T) for ^^Cu and ^^^La NQR, the spin order 
parameter measured by elastic neutron scattering, and the ^SR data, we deduce the spatial dis- 
tribution of the spin fluctuation frequency scale F and its temperature dependence. These results 
indicate that as soon as charge dynamics slow down, the spin fluctuations begin to slow dramati- 
cally with spin stiffness ^ ~ 200 K. By extending the analysis to other hole concentrations, we 
demonstrate a qualitative difference in the spatial variation of electronic states induced by slowing 
charge dynamics above and below x = |. 

76.60, 74.25-N, 74.72 
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I. INTRODUCTION 

The electronic phase diagram of high cuprates 
involves a number of different phases and crossovers, 
including an undoped antiferromagnetic phase (Neel 
state), an insulator- metal crossover, a spin-glass phase, 
a pseudo-gap phase, and an over-doped metallic regime, 
in addition to the superconducting phase. The complex- 
ity of the phase diagram makes it difhcult to identify 
the superconducting mechanism. Some of the aformen- 
tioned phases and crossovers may be irrelevant to the 
superconducting mechanism, while others may prove to 
be critical.EI Further experimental efforts are necessary 
to elucidate the properties of each phase, as well as the 
nature of transitions and crossovers between adjacent 
phases. 

A relatively new addition to the already complex elec- 
tronic phase diagram of high cuprates is the so-called 
stripe phase near the magic hole concentration a: « |. It 
has been demonstrated by elastic neutron scattering ex- 
periments on Nd co-doped La2-a;Sra;Cu04 by Tranquada, 
Sternlieb, Axe, Nakamura, and Uchida, that doped holes 
microscopically phase segregate into hole rivers form- 
ing a charge density wave (CDW).a The hole rivers be- 
come antiphase boundaries between hole-poor segments 
in which Cu spins form a short range spin density wave 
(SDW) order. Recent measurements based on neutron 
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scattering,ri NQEjl^^ ESR, 

charge tpanseort,E3cfl photo emission spectroscopyE£lE£ 
and /iSREJeS detect experimental signatures of the stripe 
instabilites in the underdoped regime. 

A major difficulty in applying NMR techniques to in- 
vestigate the stripe phase of high Tc cuprates is in the 
glassy nature of the slowing doW|H|-Qf stripes that is ob- 
served for most-jStriped cupratedaiHa (one exception be- 
ing La2Cu04+il3t3). Unhke ordinary second order phase 
transitions, the slowing of stripe fluctuations is very grad- 
ual and inhomogeneous. This means that the experi- 
mental signature of stripe freezing appears at different 
temperatures depending on the characteristic frequency 
of the experimental probe. If one uses a probe with a 
lower frequency scale, the stripe anomaly takes place at 
a lower temperature. Furthermore, the glassy freezing of 
the stripes do not show typical signatures of phase tran- 
sitions observed by NMR, such as a divergence of NMR 
spin-lattice and spin-spin relaxation rates 1 /Ti and 1 /T2 
at the on-set of freezing. That explains why earlier NMR 
studies of underdoped high Tc cuprates carried out be- 
fore the discovery of stripes by Tranquada et al. by neu- 
tron scattering failed to capture the various signatures of 
stripe instabilities in a convincing manner, even though 
there are some results thatj-seem, retrospectively, consis- 
tent with stripe anomalies. E3e3 

Very recently, we demonstrated that slowing of stripes 
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in La2_a;Sra;Cu04 with and without Nd or Eu rare earth 
co-doping and in Lai.875Bao.i25Cu04 can be easily cap- 
tured by measuring the anomalous-reduction (i.e. wipe- 
out) of the ''^Cu NQR intent^itj^'y 



In our earlier publications,ll31iil wc left some of the im- 
portant issues unexplored. First, if stripes keep slowing 
down with decreasing temperature, what kind of spin 
and charge density modulation does the stripe phase ex- 
hibit in the ground state at NMR/NQR time scales! Sec- 
ond, what is the spatial distribution of the stripe fluctu- 
ation frequency? Third, how does the hole concentra- 
tion affect the slowing of stripes? In order to address 
these issues, we now report a more detailed ^^Cu and 
"^La NQR study of Lai.875Bao.i25Cu04 and (Nd, Eu) 
co-doped La2-2;Sr2;Cu04 near the magic hole concentra- 



tion a; « g down to 350 mK. In section [II, we report 
Zeeman preturbed ^■^Cu and ^■^^La NQR lineshapes in- 
cluding detailed numerical simulations. Our preliminary 
lineshape measurements on ^'^Cu isotope |eiiriched sam- 
ples conducted for our earlier publicatiorJij did not re- 
veal any noticeable structures down to 1.7 K other than 
a single peak, in agreement with the earlier finding hjz. 
Tou etal. reported for a ^'^'^^Cu mixed isotope sample.E£l 
This made us suspect that stripes are still fluctuating 
quickly enouah. relative to NQR time scales to motion- 
ally narroi£3c3 the lineshapes at 1.7 K. Accordingly, we 
developed a top-loading '^He NMR system to conduct 
NQR lineshape measurements in a broad frequency range 
down to 350 mK in the hope that any stripe excitations 
with energy AE/ks ^1 K or higher would be suppressed. 
However, our new results at 350 mK demonstrate that 
a periodic spatial modulation of spin and charge den- 
sity hardly exists at NQR time scales even at 350 mK. 
This is in remarkable contrast with elastic neutron, x- 
ray and /xSR measurements that sucessfully detect short 
range spin and charge order at time scales that are three 
to seven orders of magnitude faster than NQR measure- 
ments. Our NQR results suggest that the short range 
stripe order is highly disordered or that extremely slow 
dynamic fluctua tion s persist even at 350 mK, or probably 
both. In section IV, we analyze the distribution of ^^^La 
nuclear relaxation ^^^1/Ti and demonstrate that below 
charge ordering, glassy slowing of the Cu spins fluctua- 
tions sets in and results in large spatial distribution of 
We also analyze the wipeout of ^^Cu and ^^^La. 
We employ the renormalized classical scaling in the non- 
linear-cr-model to fit the entire set of existing data with 
small spin stiffness < 200 K) which is reduced 

due to slowed charge dynamics. This analysis accounts 
for the temperature dependences of not only the ^^Cu 
NQR wipeout, but also the ^^^La NQR wipeout, the or- 
der parameter of spin stripes measured by elastic neutron 
scattering, the ^SR asymmetry, and the recovery of both 
^^Cu and ^^^La NQR intensity, all within a single frame- 
work. From the analysis, we estimate the spatial distribu- 
tion of stripe fluctuations as a function of temperature. 
We find that below Tcharge the glassy slowing down of 



stripes results in three orders of magnitude distribution 
in fluctuation frequencies in agreement with our analysis 
of ^■^^1/Ti. This is in contrast to the well known neutron 
results in the normal metallic state of La2_a;Sr3;Cu04, 
where the spin fluctuation energy scale hT is well-defined 
by a single value on the order of ksT, and that hV varies 
roughly in ppsportion to T as origiiially discovered by 
Keimer etal£ii (called oj /T -scalingjL^LJ or quantum crit- 
ical behaviortB). In section V, we extend the analysis to 
deduce a diagram of temperature versus stripe fluctua- 
tion frequency for various hole concentrations. 



II. EXPERIMENTAL 

All of the single phase, polycrystalline samples of 
La2-2;-y(Eu,Nd)y(Ba,Sr)2,Cu04 used in this study were 
prepared using conventional solid state procedures. We 
mix predried La203 (99.99%), EU2O3 (99.99%), Nd203 
(99.99%), BaCOa (99.95%), SrCOg (99.99%), and CuO 
(99.995%) with correct nominal compositions, and grind 
by hand with an agate mortar and pestle until an in- 
timate mixture is obtained. To make complicated line- 
shape analysis feasible, some samples were enriched with 
either ^^Cu or ^^Cu isotopes. This first grinding usually 
consumes 40-60 minutes for the 1-4 g of powder that we 
prepare. A prereaction is carried out for 20 h in a box 
furnace at 850°C followed by repeated regrindings and 
sinterings (also 20 h) at temperatures between 950°C and 
1000°C. Finally the samples are pelletized and high tem- 
perature annealed in flowing O2 gas at 1100°C to 1150°C 
for 24 h to 48 h before a slow and controlled cooling cycle 
that includes low temperature annealing at 800° C (24 h) 
and 500° C (24 h). By using a large number of grind- 
ings (typically 5-8), we acheive homogenous, high quality 
polycrystaline samples. The powder x-ray spectra reveal 
that all samples are single phased. The measured lat- 
tice parameters show systematic variation through the 
sample families as a function of hole concentration, re- 
vealing the room temperature orthorhombic to tetrago- 
nal structural transition, apd are well within the range of 
previously published data.c3 To further characterize the 
prepared samples, magnetization measurements are per- 
formed with a SQUID magnetometer at a field of 10 Oe. 
All superconducting samples show a clear diamagnetic 
response, and Tc and vokime fraction are consistent with 
previous measurements. eZI Nonsuperconducting samples 
show no Meissner shielding down to 5 K. 

All nuclear resonance data are taken with a homemade, 
phase coherent pulsed NMR spectrometer using the stan- 
dard 7r/2 — r — TT pulse sequence with a typical pulse sep- 
aration time r = 12 /xs. For measurements at and above 
1.7 K we use a commercial cryostat and measure temper- 
ature with calibrated Cernox and carbon-glass resistors. 
In order to conduct Zeeman perturbed NQR measure- 
ments in zero applied field down to 350 mK, we devel- 
oped a top-loading NMR system utilizing a commercial 
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single-shot ^He refrigerator. In those cases, temperature 
is measured via a cahbrated ruthenium oxide tempera- 
ture sensor. 

III. ZEEMAN PERTURBED NQR LINESHAPE 

A. Principles of Zeeman perturbed NQR 

In the Neel ordered state of the undoped parent com- 
pounds of high Tc cuprates, Zeeman perturbed NQR has 
been sucessfully utihzed to measure the sublattice mag- 
netization {Si) and its temperature dependence. In anti- 
ferromagnets below the Neel temperature T/v, the nuclear 
resonance frequency fi of the j-th site depends on {Hi), 
the static component of hyperfine magnetic field at that 
site. Typically, {Hi) originates from the hyperfine cou- 
pling A" with the ordered magnetic moment {Si) at the 
same site and from the transferred hyperfine coupling A^^ 
with the ordered magnetic moment {Sj) at the nearest 
neighbor sites. Since the hyperfine coupling constants 
can be determined based on NMR measurements in the 
paramagnetic state above T/v, one can utilize the nuclear 
resonance frequency fi to measure the hyperfine mag- 
netio-^ieljd {Hi), and hence the sublattice magnetization 
(S'i).E3cj This technique has been employed successfully 
in the Neel state o£fie«eral undoped parent compounds of 
high Tc cupratesErEj In tie case of undoped La2Cu04 
(Ttv = 325 K), Tsuda et alS^ find three peaks correspond- 
ing to the transitions \Iz — -f-|) '-^ \ + \),\ + \) ^ \ — \), 
and I — 5) ^ I ~ §) foi' 6ach of the ^^Cu and ^^Cu isotopes 
(both have nuclear spin / = |). Defining = A" as the 
x-component of the on-site hyperfine coupling constant 
and B = A^^ as the isotropic supertransferred hyper- 
fine coupling constant from the four nearest neighbor Cu 
siteSjtfl the hyperfine magnetic field at the nuclear site i 
can be written as 

(i/.) = A, (5.) +5^(5,) (1) 

where the sum is taken over the four nearest neighbors. 
We also use the facL-that the ordered Cu moments lie 
within the ab-plane.E3 For the present case of a two di- 
mensional square-lattice antiferromagnet, {Si) ~ ~{Sj) 
for nearest neighbors i and j. Then Eq. ^ simplifies to 

{H,) = {A^-4B){S,). (2) 

From the analysis of the three transitions of ^"^Cu zero 
field NMR (see Fig. |l|), Tsuda etal. find the internal field, 
at the Cu nucleus to be {H,) = 8.2 T in La2Cu04.E2l 
On the other hand, the hyperfine coupling constants-ajeej 
known to be A^ = 38 kOe/^s and B = 42 kOe//ii3.Ea'y 
Together, these values imply an effective moment of 
— ~ 0.63 /iB, in accordance with previous 

results of muon spin precession pijsasurements, which in- 
dicate a moment of ~0.58 /is .Ell Additionally, the fre- 
quency splitting between two satellite transitions, \Iz = 



+i) ^ 1 + 5) ^^"^ I ~ i) ^ I ~ f) ' provides information re- 
garding the nuclear quadrupole interaction tensor, which 
reflects the second derivatisje^f the electrostatic Coulomb 
potential at the Cu sites.EJil3 The splitting data may be 
used, in principle, to obtain informajtipn regarding the 
charge state of the observed nuclei Ej'EJ 

If well defined static charge and magnetic order exist 
in the striped materials, we would expect similar, highly 
structured resonance lineshapes arising from each unique 
Cu site, from which we could gain insight into the spin 
state from the measured hyperfine magnetic field Hi and 
the charge state from the nuclear quadrupole interaction 
tensor vq. Contrary to these expectations, we observed 
no evidence for well-defined periodic spin and charge den- 
sity modulations even at 350 mK, as will be described in 
the following sections. 




Hyperfine Field [T] 

FIG. 1. Plot of '^^Cu NQR (/ = |) transition frequen- 
cies as a function of hyperfine field in the ab-plane with 
an axially symmetric quadrupole tensor Vq — 34.5 MHz. 
This chart is applicable for the case of ®^Cu resonance in 
Lai.g75Bao.i25Cu04. The dashed lines mark transitions that 
occur with very low probability. The intersects of three solid 
curves and the dotted vertical line a.t H = 8.2 T correspond 
to three zerp-|field Cu NMR lines observed for La2Cu04 by 
Tsuda etalB (To be rigorous, i/^ = 31.9 MHz at 1.3 K in 
La2Cu04.) 

B. "^Cu NQR/NMR spectra in Lai.875Bao.i25Cu04 

In 1991, magnetic order was discovered in a a; = | 
doped, sample of La2-£cBaa;Cu04 using /iSR by Luke 
etal2Ji The next year, Tou etal. reported anomalous 
broadening of the low temperature ^^'^^Cu NQR spec- 
tra, consistent with Zeeman perturbed NQR.E3 Here we 
provide a comprehensive understanding of the low tem- 
perature Cu resonance lineshape in light of the progress 
that has been made in the field oL stripe physics. As 
demonstrated earlier by Hunt e<aZ£3 and Singer etal.fA 
most of the NQR response of Ba, Nd, and Eu doped ma- 
terials is semi-quantitatively identical for a; « |. In the 
following, we focus our discussion of the NQR lineshapes 
on the Lai.875Bao.i25Cu04 sample. The Nd doped mate- 
rials are ill suited for Zeeman perturbed NQR measure- 
ments at low temperatures due to the large Nd'^^ local 
moments that make the relaxation rates too fast, and the 
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Eu doped samp les show extra NQR peaks as discussed 
in section [II F 
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FIG. 2. ''^Cu NQR lineshapes of a '^^Cu isotope enriched 
Lai.875Bao.i25Cu04 sample at selected temperatures, includ- 
ing pure NQR from 290K down to 70 K, partially "wiped out" 
NQR at 55 K, 45 K and 35 K and Zeeman perturbed NQR at 
1.7 K and 350 mK. '''^Cu resonance continues below 25 MHz 
in the Zeeman pertubed case, but is obscured by the ^^^La 
NQR = ±|) |±|) transition peaked near 18 MHz which 
is more than an order of magnitude stronger in intensity (not 
shown). Note the scale multipliers when comparing magni- 
tudes. The solid curves are guides for the eye. The 1.7 K and 
350 mK curves are identical except for overall magnitude. 

Selected ^^Cu NQR lineshapes for ^^Cu isotope en- 
riched Lai.875Bao.i25Cu04 are shown in Fig. |[ At 
high temperatiires we observe the well known A and 
B NQR lines.E3 The relative intensity of the B line is 
the same as the doped concentration x of the impurity 
ions Ba^"*" and Sc^r-ii over a broad concentration range 
0.02 < X < O.SoMeI Furthermore, the B hne is known 
to be observable even at elevated temperatures as high as 
800K.l3 Accordingly, it is unlikely that the B line origi- 
nates from electronic effects such as phase separation or 
stripes. The most likely scenario is that the B line arises 
from Cu sites directly across the apical O from the doped 
Ba position and the lattice contribution to the electric 
field gradient is different due to the presence of a Ba^+ 
ion rather than La^+. Needless to say, this line assign- 
ment of the B line does not exclude the possibility that 
the B site has a somewhat different local electronic state 



as compared to the surrounding A sites. In fact 
is known to be somewhat smaller for the B site as shown 
by Yoshimura et. al.E3 In Fig. ^ we also show that Eu 
codoping results in a new structure at the low frequency 
side of the A line. We attribute this new structure to the 
Cu site sitting across the apical O from the Eu'^+ ions. 
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FIG. 3. Integrated '^^Cu NQR intensity above 25 MHz 
for Lai.875Bao.i25Cu04 (♦), Lai.68Euo.2oSro.i2Cu04 (a) and 
Lai.4gNdo.4oSro.i2Cu04 (•) scaled for the Boltzman factor 
and T2 spin-echo decay. In addition, we plot the charge or- 
der parameter (o) deiduced from neutron scattering on the 
Nd codoped material.B Within the broadened hyperfine field 
model in Fig. ^f, the maximum recovery of signal above 
25 MHz is expected to be 61%, as indicated by the arrow. 
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FIG. 4. The quadrupole resonance frequency, ^'"^uq, ob- 
served at fixed pulse separation time r = 12 fisec for the A 
line (o) and B line (•) of a ^'^Cu isotope enriched sample of 
Lai.875Bao.i25Cu04, as well as the peak frequency of the low 
temperauture, Zeeman perturbed NQR spectra (A). 

The temperature dependence of the integrated ^^Cu 
NQR intensity is presented in Fig. ^. Since the occupa- 
tion probability of nuclear spin levels is determined by 
the Boltzman factor, the bare NQR intensity is inversely 
proportional to temperature T. Accordingly, we multiply 
the bare NQR intensity by T to correct for the Boltzman 
factor. We also calibrate the intensity of the observable 
NQR signal by extrapolating the observed spin echo de- 
cay to 2r = 0. The NQR intensity thus deduced corre- 
sponds to the number of nuclear spins whose resonance 
conditions are well defined and whose relaxation times 
^^Ti and ^^T2 are long enough to be observed experimen- 
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tally within the time domain of pulsed NQR measure- 
ments, T]\fQji ~ 20 /isec. Missing signal intensity arises 
because at ileast one of these conditions is violated for the 
lost signal.E3'EJ We note that these are standard proce- 
dures that have also been applied to the NMR investiga- 
tion of the SDW state in electron doped Nd2-xCej;Cu04 
by Kambe et. al.xl Between room temperature and 70 K 
the lines retain constant intensity, as expected, because 
the number of ^^Cu nuclear spins detected by NQR mea- 
surements is unchanged. The peak frequency, vq{A), 
displays only a small temperature dependence down to 
Tcharge ~ 70 K as shown in Fig. IJ. The weak tempera- 
ture dependence of vq (A) down to 70 K can be accounted, 
for entirely based on slight tilting of CuOe octahedra.tll 
However, as the temperature is decreased below 70 K, 
the temperature dependence of both i^q (A) and vq (B) 
changes, and both A and B lines display an identical 
drop in intensity and both nearly disappear by 15 K.til 
This drop in signal intensity is "wipeout" . As previously 
noted by Singer etai, the estimation of the NQR inten- 
sity close to the onset of wipeout is made more difficult 
due to the change in the form of the spin-echo decay. 
Given the uncertainties in fitting the spin echo decay, 
it is ininossible to define T^qr with accuracy exceeding 
±10%.llll Within these uncertainties, the onset of wipe- 
out, Tnqb., is indistinguishable from Tcharge for x ~ |. 
The wipeout and subsequent recovery of the intensity be- 
low 15 K will be analyzed in detail in section IV. As the 
temperature is decreased below 8 K, the resonant signal 
begins to reemerge. This new line is shown in the bottom 
panel of Fig. ^. We note that the character of the line is 
entirely altered. The peak has shifted down in frequency 
from 35.5 MHz to 31 MHz. Additionally, the full width 
at half maximum, which slowly increases with decreasing 
temperature from 2.0 MHz at 300K to 2.6 MHz at 70 K to 
3.2 MHz at 15 K, shows sudden and dramatic broadening 
to more than 12 MHz at 350 mK with the emergence of a 
large tail towards higher frequencies. Note that the spec- 
tra at 1.7 K and 350 mK are of identical shape and differ 
only in magnitude. In fact, the curves through those 
data sets in Fig. ^ are identical except for vertical scal- 
ing. Based on the analysis presented in the following two 
sections, we demonstrate that these changes are consis- 
tent with a broad distribution in the hyperfine magnetic 
field (Hi) averaged over NQR time scales. 



the descriptions of the spin and charge density modula- 
tion based on this model and similar models as shown in 
Fig. ^-e, (which are consistent with the neutron data) 
are valid for the entire sample at the much slower NQR 
timescale of ^ 20 /isec, we would expect inequivalent Cu 
sites. In the case of the model of Fig. we expect there 
to be three inequivalent Cu positions: one site in the 
middle of the magnetic region, one at the edge and one 
in the charge rich river. As a first approximation, we ex- 
pect the two magnetic sites at the center and edge of the 
stripe to yield three NMR lines each, while the nonmag- 
netic boundary site should produce a single NQR line. 
Thus, we might expect a total of at least seven resonance 
lines. These lines are not identifiable in the data. 

Given that the 350 mK lineshape shown in Fig. || is 
lacking the expected structure, we turn to numerical sim- 
ulations by solving the Hamiltonian exactly. Using / as 
the nuclear spin, rj as the asymmetry parameter of the 
electric field gradient (EFG), 7„ as the gyromagnetic ra- 
tio, h as Planck's constant, and H as the static hypecfine 
magnetic field, the Hamiltonian can be written asEjO 

_ + 1) ^ l^(^2 + ^2 )| _ ^^f^j^ . 

(3) 

where the ^'^Cu nuclear spin is / = | and the quadrupole 
resonance frequency ^^i^q, is related to the electron 
charge e, the maximum of the EFG q, and the ^^Cu 
quadrupole moment ^'^Q, by 



63, 3e^g 



2hI{I+l) 



(4) 



Because of local axial symmetry of the Cu sites, the elec- 
tric field gradient is nearly axially symmetric in both 
undoped and hole doped Cu02 pUfies. Additionally, 
as originally shown by Ohsugi et aLEjc3 the ^^^La NQR 
linebroadening takes place only for the | ± i) ^ | ± |) 
transiti on in the low temperature striped phase (also see 
section [II G). This indicates that the internal field lies 
primarily in the Cu02 plane for x ~ ^ E3S Under these 
conditions, the Hamiltonian from Eq. His greatly simpli- 
fied and can be written in matrix form, where we chose 
the basis to be the eigentates of the operator Iz, 



C. Numerical simulations in standard stripe models 

At the time of the discovery of stripes in Nd doped 
Lai.48Ndo.4oSro.i2Cu04 by neutron scattering,0 Tran- 
quada et al. proposed a simple model valid at the short 
timescale of neutron scattering (1 meV w 2.4 x 10^^ Hz, 
corresponding to 4psec) which consists of sets of three 
magnetic rows (a "stripe" ) separated by a charge rich 
river which serves as an antiphase boundary for the anti- 
ferromagnetic order. This scenario is shown in Fig. If 
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(5) 

Between the four eigenstates of the Hamiltonian, there 
are six possible transitions corresponding to six resonant 
frequencies as shown in Fig. ^ Some of these transi- 
tions are nearly forbidden by the |Am| = 1 selection rule 
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and exist only due to the unusual mixing of states that 
occurs in the low and intermediate field regions where 
InH < vq. Thus we must consider the probabilifas 
for each transition to occur. Following G.H. Muha,0 
we treat the small experimentally induced R.F. field, 
Hi (^lOOG), as a perturbation and use the the stan- 
dard time-dependent expression for the transition rate 
between two states n and n' 



|(n'|I-Hi|n) 



(6) 



Since we are working on polycrystalline samples where 
the direction of Hi with respect to the crystal axis is 
random, we perform a powder average to obtain the 
probability for each transition. These probabilities cor- 
respond to the expected intensity for each resonance line 
in Fig. P Full lineshape simulations are computed nu- 
merically by summing contributions from distributions 
in both the hyperfine field and in the quadrupole fre- 
quency. If the charge density modulation is sizable, ^'^vq 
may vary depending on the site, and the axial symme- 
try of the EFG tensor can break down. However, the 
structureless broad resonance line observed below 8 K 
suggests that the spatial modulation of spin and charge is 
smooth at the NQR time scale, r^qq^ « 10~^ sec. Hence 
we believe that approximating the quadrupole tensor as 
axially symmetric is acceptable in our simulations. Ex- 
trapolating from observations above 70 K before wipeout 
begins we find ^^vq{A) = 34.5 MHz for the A site and 
63jyQ(S) = 41.8MHz for the B site. Note that, although 
the apparent peak frequency of the A line increases some- 
what below 70 K as shown in Fig. I 63j,q(A) = 35.5MHz 
at 10 K is achieved only for several percent of the sample. 
It is not initially clear what distribution of quadrupole 
frequencies to employ at low temperatures, where we ex- 
pect charge order to induce added width. In this section, 
we take both the A and B lines to be Gaussians, each 
with a width of 5.6 MHz (full width at half maximum), 
corresponding to twice the value extrapolated from the 
A Hue data above 70 K. We discuss the effects of differ- 



ing amounts of quadrupole broadening in section |IIIE . 
To allow direct comparisons to the data, the simulated 
intensities are weighted by the square of the resonant 
frequency to account for experimental sensitivity.cil 

The next step is to estimate the distribution of the hy- 
perfine magnetic field in various stripe models. First we 
examine the aforementioned model that was initially pro- 
posed by Tranquada et al. ,□ which is shown schematically 
at the left of Fig. ^a. The Cu site in the middle of the 
magnetic "stripe" is locally in the same enviornment as 
in the antiferrmagnetic Neel state, so we expect the field 
at the nucleus to be ijc'^"*^'' = \A^- AB\^i^}f 
an effective magnetic moment of /ie// = 0.3 /is 
Inserting this value and the hyperfine coupling constants 
discussed above in section HI A, we estimate = 
130kOe/MB X 0.3 = 3.9 T. For the Cu site at the 
edge of the stripe, there are only three magnetic nearest 
neighbor sites that give rise to supertransferred hyperfine 



fields, so the expected field is iJ^'^f^ = \Ax — 3B\iJ,eff = 
2.6 T. The last Cu site is within the antiphase boundary, 
where we expect oppositely oriented magnetic contribu- 
tions from two nearest neighbors which cancel, yielding 

^boundary ^ q ^^ot thcSC valuCS of i7<=enter^ jjedge 

and ffboundary ^j^g distribution of internal field at the 
top center of Fig. |^, where the vertical scale is determined 
by the relative numbers of each Cu site, which is the in- 
tegral of each peak. The widths are chosen as depicted in 
the figure, where for H > the widths are proportional to 
the hyperfine field H as would be expected for magnetic 
broadening. The numerically computed ^"^Cu lineshape 
is shown at the right of Fig. |^a. Agreement between the 
data and the simulated lineshape is poor. Specifically, all 
of the structures are too sharp, and even the largest peak 
is located at ~ 37 MHz, above the 31 MHz peak seen in 
the data. 

In addition to the original model proposed by Tran- 
quada et al. , we examine site centered and bond centered 
commensurate sinusoidal spin density waves and an in- 
commensurate sinusoidal spin density wave, as well as 
a grid model (site centered sinusoidal spin density wave 
in two directions) as shown in Fig. ^-e. In each case 
/ie// = 0.3 i^B represents the magnitude of the maximum 
moment. For each model, we follow the same procedure 
as outlined above to construct a hyperfine distribution 
and compute an expected ^'^Cu lineshape. In each of the 
simulated lineshapes there is too much structure for a 
good match to the smooth data. 

Even though agreement between the experimental 
spectra and the numerically simulated lineshapes is poor, 
we notice that, in Fig. ^ and there is substantial 
simulated signal at the experimental peak frequency of 
31 MHz. These contributions do not arise from the de- 
tails of the hyperfine field distribution, but rather from 
the nature of the transition energy levels themselves. If 
we look back to Fig. |l], we note the two transitions that 
start at vq = 34.5 MHz and move to lower frequencies 
with increasing field. These lines achieve a broad mini- 
mum at 31 MHz. Because of the high density of states 
that corresponds to resonant frequencies near 31 MHz, 
any continuous distribution of internal fields that includes 
contributions in the broad window from 0.5 T to 2.5 T 
will give rise to a lineshape with a 31 MHz peak. With 
this in mind we examine broad distributions of hyperfine 
field, and by trial and error we achieve a good fit to ex- 
perimental data as shown in Fig. ||f. In section HI E we 
will return to the issue of the low temperature lineshape 
to understand the magnetic state of the Cu-0 plane at 
the slow timescales of our pulsed NQR experiments. 
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FIG. 5. Numerical simulations of the Cu spectra based 
on several microscopic stripe models. The first column gives 
a schematic diagram of the orientation and magnitude of the 
Cu magnetic moments along the Cu-0 bond direction. Ar- 
rows and open circles represent ordered Cu moments and a 
hole-rich anti-phase boundary without ordered moment, re- 
spectively. The middle column displays the expected distri- 
bution of the hyperfine magnetic field at the Cu sites. The 
right column contains the corresponding simulated lineshape. 
(a)-(c) represent several commensurate stripe models, and 
an incommensurate spin density wave is presented in (d). 
Model (e) is a two-dimensional sinusoidal grid structure in 
the Cu02 plane and (f) depicts a simple, broad distribu- 
tion of hyperfine field that provides a good fit to the 350 mK 
Lai.875Bao.i25Cu04 data, which is included as open circles. 



D. Recovered ^^Cu NQR intensity at 350 mK 

So far we have considered only the lineshape. Another 
important point that needs to be addressed is the inte- 
grated NQR intensity at 350 mK. The accurate measure- 
ment of the recovered integrated NQR intensity below 
8K is not easy. The ^"^Cu NQR lineshape measurements 
above 10 K can be performed in one single scan over fre- 
quency without changing the experimental set up, but 
the Zeeman perturbed NQR lineshape is so broad that 
we need to replace the RF coil of the tuned LC-circuit 
in our NMR probe three or four times to cover the en- 
tire frequency range from 20 to 90 MHz. This enhances 



uncertainties in our estimate of the integrated intensity, 
and this is the primary reason why we have large error 
bars for the integrated intensity below 8 K in Fig. ||. In 
addition, the lineshape below 25 MHz is superposed by 
an order of magnitude stronger ^^^La NQR line centered 
at 18 MHz. This does not allow us to measure the ^^Cu 
Zeeman perturbed NQR intensity below 25 MHz. Above 
70 K, earlier high field measurements for an aligned pow- 
der sanmle guarantee that there is no NQR line below 
25MHzO Within the large uncertainties, the recovered 
signal intensity integrated above 25 MHz represents only 
50^20*^ of the value above 70 K. However, it is important 
to notice that for hyperfine fields ranging up to ~ 4 T, 
there are two branches of Zeeman perturbed NQR tran- 
sitions below 25 MHz in Fig. |l|. Within our model calcula- 
tions of Fig. that correctly reproduce the lineshape, we 
can estimate the fraction of integrated NQR intensity ex- 
pected to appear below 25MHz, which is 39%. It should 
be emphasized that the echo intensity for the simulated 
lineshape plotted in Fig. ||f is scaled by multiplying by the 
square of the frequency to take into account the experi- 
mental sensitivity. Accordingly, the apparent integrated 
intensity below 25 MHz incorrectly appears to be a tiny 
fraction of the whole lineshape. In other words, the max- 
imum integrated intensity expected to reemmerge above 
25 MHz is 61% of that at 70 K, as shown by a horizontal 
arrow in Fig. ^. This means that the experimentally ob- 
served intensity at 350 mK is effectively 82^33% of that 
at 70 K, accounting for nearly the full intensity from the 
sample. This rules out the possibility that the descrep- 
ancy between the experimentally observed Zeeman per- 
turbed NQR lineshape and theoretical calcualtions based 
on stripe models originates from missing intensity. 



E. Numerical simulations based on stripe models 
including broadening and motionally narrowing 

The good agreement between the data at 350 mK and 
the simulated lineshape and its integrated intensity based 
on the broad distribution of Fig. &, indicates that the 
hyperfine field averaged over NMR time scales and over 
the entire sample should extend up to about 3 T with 
substantial weight in the region below 2 T. Motivated by 
the success of the fit for the broad distribution of hyper- 
fine fields, we return our attention to the stripe models 
of Fig. ^-e to examine the possibility that the internal 
hyperfine field corresponding to the stripe models may 
be broadened, either by disorder within the sample or 
by motional averaging effects, in such a manner as to be 
consistent with the ^^Cu NQR spectra. 

In Fig. 1^ we develop the hyperfine field distribution 
expected from an incommensurate-, spin density wave, 
favored by a recent /iSR studypj including varying 
amounts of broadening, and calculate the correspond- 
ing simulated ^^Cu hneshape. We take the magnitude 
of the maximum spin in different regions of the sample 
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to be (S), which may be distributed due to static inho- 
mogenieties. We numerically simulate the hyperfine field 
based on a row of 10,000 electron spins of magnitude 
(5„) = {S) X sin(27rn/(8 + S)) with \S\ < 1, (5 ^ 0, taken 
to be aligned along the x axis and labeled by position 
n. ThjB spatial periodicity is that expected from neutron 
data.H The internal field at each site is computed with 
the use of Eq. |l|, giving = A^{Sn) - 2B{-Sn) - 

B{Sn+i) — B{Sn-i) where we assume Neel order in the y 
direction and use the values of and B as given above. 
Distributions arising from different segments of the sam- 
ple (different {S) values) are combined as per the distri- 
bution in (5*). In Fig. ^-c we show three simulations 
based on distributions in the magnetic moment that are 
centered at 0.3 /^b. With increasing width in {S), both 
the hyperfine distributions and the simulated lineshapes 
become more smooth. However, the smoothed lineshaped 
are shifted to too high a frequency to match the experi- 
mental data. Thus we rule out the possibility of a static 
incommensurate spin density wave at NQR time scales 
with maximum magnetic moment of 0.3 /i^, even with 
possible broadening of the value of the moment. 
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FIG. 6. Models of incommensurate spin density waves 
based on a distribution of maximal moment in different seg- 
ments of the sample as given in the leftmost column. Corre- 
sponding to these distributions we plot the hyperfine field at 
the nuclear site in the middle column and the simulated ^''Cu 
lineshape at the right. In models (a-c) and (e) we assume 
quadrupole width of 5.6 MHz; in (d) that width is taken to 
be 7.5 MHz and in (f) 3.7 MHz. The 350 mK spectra data is 
shown as open circles. 

Next we consider the possibilities of a smaller average 
ordered moment and of a different amount of quadrupo- 
lar broadening. The simulations of Fig. ^-f represent 
good fits to the data based on average ordered mo- 
ments of 0.10, 0.15 and 0.20 paired with quadrupole 
widths of 7.5, 5.6 and 3.7 MHz respectively. Without 
this extra quadrupole broadening the contributions of 
the A and B lines give rise to distinct peaks in the sim- 
ulated lineshape, in contradiction with our experimental 
findings. In the homogeniously doped Cu02 planes of 
La2-2;Ba2:Cu04 with x > the doping dependence of 



vq {A) is known.Ea If we imagine the system to consist of 
many patches, each of which being homgeniously doped, 
the extra width in i/q (A) beyond the extrapolated value 
from T > Tcharge (2.8 MHz) would correspond to dis- 
tributions in X of approximate width Aa; ~ 0.30, 0.18 
and 0.06 respectively. There is a trade off between the 
broadening caused by magnetic effects and that caused 
by inherent distribution of the quadrupole resonance. 
However, the simulated lineshape of Fig. |6^ best fits 
the data, indicating that within that incommensurate 
spin density wave model a distribution of (S) centered 
at 0.15 fiB coupled with moderate quadrupolar broaden- 
ing corresponding to Ax ~ 0.18 provides the best fit to 
the data, although there is a range of parameters ((5) 
and quadrupole width) that provide a reasonable fit to 
the data, as indicated in Fig. p|d,f. Of course this model 
assumes patches of homogenious doping rather than mi- 
croscopic structure (stripes) which represents an alterna- 
tive, and perhaps more likely, path to increased width in 

How can we-jreconcile this result with the conclusion of 
Kojima etalE3 based on their /iSR data? Kojima etal. 
observed well defined precession of positive muons con- 
sistent with an incommensurate spin density wajv« with 
well defined maximum ordered moment of 0.3yLtB£2l If we 
recall that muon spin precession persists for only about 
0.5 /isec, while the NQR spectrum is obtained by spin 
echo measurements that detect hyperfine field averaged 
over 20 /isec, one possible explanation is that the rela- 
tively large and well-d efined m oment of 0.3 /is observed 
by /iSR measurement£3'ESEil is still slowly fluctuating 
during the duration of a Zeeman perturbed NQR mea- 
surent. Since the hyperfine field is a vector quantity, the 
averaging process will always reduce the apparent, av- 
eraged field. Thus, the effect of motional averaging will 
be to pull all instantaneous values of the hyperfine field 
towards zero, reducing the effective values. This may be 
the key to resolving the apparant discrepancy between 
our NQR findings and the /iSR measurements. We rp^^ 
that resistivity never diverges in the stripe phaseJHU'EHI 
hence we believe that mobile holes disrupt spin order. 
We note that the motional averaging scenario is sup- 
ported by the large value of at 350 mK (see 
Fig. which indicates that low frequency spin fluc- 
tuations persist even at our experimental base temper- 
ature. We also note that the motional averaging effects 
in the present case of Zeeman perturbed, JMjQR is some- 
what different from the textbook caseaidS of motional 
narrowing and exchange narrowing. In those cases, one 
applies a strong and static external magnetic field to ob- 
serve NMR signals. The applied field defines the Larmor 
precession frequency of nuclear spins, and the dipole line 
brodening is only a small perturbation. In the present 
case, however, the fluctuating hyperfine field caused by 
freezing Cu moments is comparable to or even greater 
than the unperturbed NQR frequency. Thus, as long 
as very slow fluctuations exist in the hyperfine field, the 
Larmor precession itself is not well defined. Accordingly 
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we expect that our Zeeman perturbed NQR spectra is 
extremely sensitive to very slow fluctuations of the hy- 
perfine field. It is worth noting that our case resembles 
the motional narrowing conside|:ed byKubo and Toyabe 
for zero or small field resonanceE3 that has beeo applied 
to the analysis of /iSR data in glassy systemsj£3 

Although we have focused on the incommensurate spin 
density wave model here, the possibilities of static broad- 
ening and motional averaging can be applied to all of the 
models presented in Fig. In each case the sharp struc- 
tures can be smoothed out by inhomogeneities within the 
sample or by motional effects, and the hyperfine field dis- 
tributions can be made to be similar to that of Fig. |^, 
and hence with the experimental data. Thus, we cannot 
exclude any of the stripe models in Fig. ^ based on our 
simulations. 



F. ^^Cu and ^^Cu resonance spectra in 
Lai.68Euo.2oSro.i2Cu04 and Lai.64Euo.2oSro.i6Cu04 

In Fig. we show the Zeeman perturbed NQR spectum 
for two samples of Lai.68Euo.2oSro.i2Cu04, one enriched 
with the ^^Cu isotope and the other with ^^Cu. The 
peak frequency of the major peak is 34±0.75MHz for 
the ^^Cu enriched sample and 31±0.75MHz for the sam- 
ple enriched with ^^Cu. The peak frequency of the ^^Cu 
sample is 3±1.5MHz higher than the 31MHz observed 
for ^^Cu in Lai.875Bao.i25Cu04. The shift of the Zeeman 
perturbed NQR peak between these two materials arises 
primarily from the difference in the intrinsic quadrupole 
frequencies, which is 2.0 MHz at 70 K. Additionally, the 
greater proximity of the B line in the Eu codoped mate- 
rial (see the inset of Fig. |^) contributes to a higher peak 
frequency in the Zeeman perturbed spectra, whereas the 
greater separation of the B line in Lai.875Bao.i25Cu04 
results in increased width. 

The only extra feature visible in these lineshapes 
for the Eu doped materials as compared with those of 
Lai.875Bao.i25Cu04 (recall Fig. ||), is the peak at 62 MHz. 
We also observed extremely weak signals that continue 
above 80 MHz. We believe that these signals arise from 
NQR of the Eu nuclei themselves. We note that the peak 
at 62 MHz is observed only in the Eu doped materials and 
is independent of the Cu isotope present in the ^^'^^Cu 
isotope enriched samples. The relaxation time at the 
62 MHz peak is also slow, allowing us to use increased 
delay time r to separate it as shown by open triangles in 
Fig. 0a. The reduced relaxation rate is consistent with 
the small hyperfine field expected at the Eu site, which 
is substituted at the La position. Furthermore, if we as- 
sume that the electric field gradient is the same at the 
Eu and La positions, use the ratio of the quadrupole mo- 
ments for the two nuclei {^^^Q/^^^Q = 4.5), and employ 
Eq. m we find that the expected peak position of ^^^Eu 
NQR is ^^^i^Q — 59.5 MHz, very similar to our observa- 
tion. Finally, we note that the -'^^^Eu peak is still observed 



at 77 K, strongly indicating that there is no connection 
with the low temperature Cu spectra. 
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FIG. 7. Resonance spectra at 1.7 K with a pulse separa- 
tion time of r = 12 /is for single isotope enriched samples of 
Lai.8-:.Euo.2oSr^Cu04 (''^Cu, a; = 0.12, o; '^^Cu, j; = 0.12, •; 
^^Cu, a; = 0.16, x). Data taken with longer r = 100/is on the 
"^Cu enriched sample is marked by A. The solid curves rep- 
resent numerical simulations of the x = 0.12 data as described 
in the text, while the dotted curve is a guide for the eye. In- 
set: 100 K lineshapes for the '''^Cu enriched samples using the 
same symbols as above. In the inset the solid curves repre- 
sent the sum of three Gaussians representing the A, B and C 
peaks, which are individually depicted by dashed curves for 
the X = 0.12 sample. For both cases the intensity ratio A:B;C 
is consistent with the expected form 0.8-a;:x:0.2 from Cu 
sites away from all dopants, closest to the Sr ions and next to 
the Eu dopants, respectively. 

Disregarding the Eu NQR signal, we focus on the 
Cu resonance in Lai.gsEuQ 2oSro 12CUO4 below 55 MHz. 
We repeat the same lineshape analysis as discussed 
above in section p^II C| for the Eu doped materials. We 
substitute ^^vgiA) = 36.2MHz, ^^uq{B) = 40.2MHz 
and ^^i^q(C) = 34.5 MHz as extrapolated from tem- 
peratures above the onset of wipeout and include the 
slightly greater (~ 25% greater than the 5.6 MHz used 
in Fig. 1^) width in quadrupole frequencies that is char- 
acteristic of the Eu doped materials. We employ the 
same distribution of hyperfine magnetic field as above 
for Lai.875Bao.i25Cu04 (see Fig. ||e). For the ^^Cu iso- 
tope enriched sample we use the ratio of qudrupole mo- 
ments i^^Q/^^Q = 0.924) and the gyromagnetic ra- 
tios (^^7/®'^7ii — 1.071) between isotopes to determine 
^^i'q{A,B) and to scale the response to the hyperfine 
field. In this way, there are no free parameters in the fit 
of either the ^^Cu or the ^^Cu spectra. The numerically 
calculated lineshapes are shown in Fig. ^ as solid lines. 
The good agreement between these simulated lineshapes 
with the experimental data for both isotopes supports the 
Cu spectra analysis that has been presented throughout 
section 



HI 



In order to investigate the influence of different hole 
concentrations on Zeeman perturbed ^^Cu NQR spec- 
tra, we carried out preliminary measurements in ^^Cu 
isotope enriched Lai.64Euo.2oSro.i6Cu04. The recovered 
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integrated intensity at 1.7 K is more than a factor of two 
smaller than in Lai.68Euo.2oSro.i2C«£)4 (See Fig. [lo|c), 
in agreement with Teitelbaum etalS3 The overall line- 
shape presented in Fig. is very similar to the case of 
Lai.68Euo.2oSro.i2Cu04, but the high frequency tail is 
smaller and the width of the main peak 10.5±0.5MHz, is 
narrower than 12.5±0.5MHz in Lai.68Euo.2oSro.i2Cu04. 
Within the lineshape analysis presented above in section 
III, this change in width implies that the maximum mag- 
nitude of the hyperfine field in Lai.64Euo.2oSro.i6Cu04 is 
5-10% smaller than in Lai.68Euo.2oSro.i2Cu04 at 1.7 K. 
There are a few possible explanations for the difference 
in hyperfine field. First, siijec-ateipes are less stable at 
hole concentrations x ^ i ^tl there may be ad- 
ditional fluctuations for those concentrations that may 
motionally narrow the main peak. This picture is con- 
sistent with the fact that at 1.7 K, a smaller fraction of 
the Cu signal intensity reemerges as Zeeman perturbed 
NQR in the samples away from a; « Second, it is 
possible that the increased disorder caused by the ad- 
ditional holes may increasingly frustrate the spins, and 
reduce the magnitude of ordered moments. We also note 
that the peak is slightly higher in the cc = 0.16 material 
as compared to x = 0.12, because vq{A) increases with 
hole doping. 

Before closing this section we would like to com- 
ment breifly on two recent reports of low tempera- 
ture ^^Cu NQR spc&U-a in Lai.8-2;Euo.2oSr2:Cu04llj and 
La1.875Bao.125CuO4.l23 Both of these publications show 
broad lineshapes superposed with what appears to be 
NQR signal from the high temperature A and B lines 
that have survived to low temperature. We have seen 
comparable results on a poorly annealed sample of 
Lai.875Bao.i25Cu04 in which wipeout was not complete. 
After careful reannealing of this sample, wipeout became 
complete and the A and B lines disappeared from the 
low temperature spectra. We note that potential prob- 
lems with the lineshapes does not diminish the primary 
finding of Teitel'baum et al. regarding the doping depen- 
dence of the recovered signal, i.e. that the intensity of 
the recovered signal at 1.3 K in Lai.8-2:Euo.2oSrj:Cu04 is 
greatest for x ~ ^ indicating the increased stability of 
that doping level. Ej This conclusion is in agreement with 
our earlier finding that the temperature dependence of 
wipeoU|t-effects in Lai.6-a:Ndo.4oSra;Cu04 is sharpest for 
x « i.EJ It is also consistent with another earlier find- 

o 

ing based on hyperfine broadening of the ^^^La NQR line 
at 1.3 K by Ohsugi etal. indicating that stripes are most 
stable near x ~ 0.115 in La2_a;Sr2,Cu04.E3 



G. ^^"La NQR lineshapes in Lai 875Bao.l25 CUO4 

To this point we have discussed only the results of 
63S5(j^ NQR, but resonance studies-of the ^^^La site can 
also provide additional insight .LjO'LJ Since the ^'^^La nu- 
cleus has spin / = |, the ^^^La NQR spectra of a sin- 



gle site consists of three peaks, roughly corresponding 
to 1, 2 and 3 times ^^^i'q- With the addition of a per- 
turbing magnetic field, the three resonance lines are split 
in a manner that depends on the details of the direc- 
tion and the magnitude of the hyperfine field.Ea In ear- 
lier ^^^La NQR studies of the low temperature striped 
phase of La2_2;Sr2;Cu04, Ohsugi etal. found that the 
I ± i) <-> I ± |) transitions at ~6MHz broadens for 
X « i at 1.3 K, while the | ± |) ^ I ^ |) transitiDna 
at ^ 18 MHz do not show magnetic line broadeningoa 
This implies that the tinjp-averaged hyperfine field lies 
within the Cu02 plane.EdS In this case, the extra line 
brodening is roughly equal to A-jnHab, where Hat is the 
hyperfine field. 
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FIG. 



^^^La linewidth as a function of temperature for 
the I ± |) ^ I ± §) transition at a frequency ^^^lyg « 6 MHz 
for Lai.875Bao.i25Cu04 (•) and Lai.68Euo.2oSro.i2Cu04 (o). 



Our goal is to use the broadening that we observe in 
the I ± i) ^ I ± |> ^^^La NQR transitio n as a test 
of the models presented above in sections [II C . We 



present the temperature dependence of the linewidth of 



the I ± 



±1) 



139t 



for 



2/ - . 1-^2/ NQR transition in Fig, 

Lai.875Bao.i25Cu04. Lineshape data from the same sam- 
ple is presented in Fig. ||. The linewidth broadens from 
l.OMHz at 20K to 1.9MHz at 1.7K, but does not dis- 
play the clea&^pJitting that is observed in the Neel state 



of La2Cu04 



I Similar temperature dependence was 



found in Lai.68Euo.2oSro.i2Cu04 as shown in Fig. ^. In- 
crease of linewidth at a somewhat higher temperature 
implies that slowing of spin fluctuations is more rapid 
in Lai.68Euo.2oSro.i2Cu04 than in Lai.875Bao.i25Cu04, 
consistent with the larger ^^^l/TiT at the same temper- 
ature in the former material, as seen in Fig. |l^. In a 
recent work, Teitel'baum et. al. also found-similar behav- 
ior for a Nd doped sample with x — 0.12.Ea 

To test the consistency of the ^■^^La line broadening 
with the models proposed above to account for the ^^Cu 
Zeeman perturbed lineshape, we un dertake an analysis 
very similar to that of sections III C . Again we perform 
an exact diagonalization of the Hamiltonian (Eq. |^), but 
this time for spin I = employing the same value of 
J] and the same small tilting of spin orientation as de- 
duced by Nishihara etal.E! We confirmed that our con- 
clusions do not depend on small variations in these ad- 
ditional parameters. In the right hand column of Fig. ||, 
we present three model hyperfine field distributions in- 
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eluding (a) and (b) which are the same hyperfine field 
distributions as used in the ^^Cu Hneshape analysis as 
shown in Fig. |^a and ||f. To avoid confusion, the mag- 
nitude of the hyperfine field is still specified by its value 
at the ^■^Cu site, although its value is greatly decreased 
at the ^'^^La position because of the two orders of mag- 
nitude smaller hyperfine coupling constant. Specifically, 
we take the ratio of the hyperfine field at the ^'^^La and 
^^Cu sites to be the samaas in undoped La2Cu04, which 
is 0.1T/8.2T=0.012MEa This means that if there is 4T 
of hyperfine field at a Cu site, there is 0.048 T of hyper- 
fine field at nearby La sites. We must caution that this 
ratio can, in principle, be dependent on the exact mate- 
rial chosen and especially on low temperature structural 
transitions which alter the mixing of atomic orbitals. In 
addition, if there is a strong modulation of the hyper- 
fine field at nearby Cu sites, the La ions may experience 
an average of those values. However, we find that the 
use of this ratio 0.012 adequately describes our data. Fi- 
nally, we take the 20 K "9La NQR data of Fig. | and 
numerically simulate the added width that arises from 
the static hyperfine field distributions given by the dis- 
tributions (a-c) of Fig. ^. These simulated low temper- 
ature lineshapes are depicted as solid lines in the main 
panel of Fig ^. We note that the differences between 
the simulated lineshapes corresponding to the different 
hyperfine field distributions are small, indicating that it 
is impossible to choose one over another based on the 
^^^La NQR measurement alone. However, we also note 
that all are consistent with the observed broadening of 
the ^•^^La NQR lineshape, giving added weight to the 
correctness of our explanation of the ^^'^^Cu and ^^^La 
spectra in Lai.875Bao.i25Cu04 as presented throughout 
this section. 
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^^^La NQR spectra of Lai.875Bao.i25Cu04 for the 
± |) transition at 20 K (o) and 1.7K (•) in the 



FIG 

main plot to the left. The solid lines represent simulated line- 
shapes that arise from three different distributions of internal 
hyperfine magnetic field as discussed in the text. The hyper- 
fine distributions are depicted in the small plots to the right, 
and are labeled (a-c) to match the corresponding simulated 
lineshape. 



H. Summary of Zeeman perturbed NQR lineshape 
measurements and analysis 

We have demonstrated that for Lai.875Bao.i25Cu04, 
Lai.68Euo.2oSro.i2Cu04 and Lai.64Euo.2oSro.i6Cu04, the 
^^Cu NQR signal that is wiped out at high temepratures 
reemerges below ~ 8 K as a single, broad peak with a 
lower peak frequency. The broad linewidth is caused by 
hyperfine magnetic fields from frozen Cu magnetic mo- 
ments that are nearly static over NQR time scales. The 
hyperfine field H is primarily within the Cu02 planes 
and the hyperfine fields extend up to 3 — 4 T. The upper- 
bound of the field distribution corresponds to the maxi- 
mum frozen Cu magnetic moments of 0.2/iB — O.S/is, in 
rough agreement with estimates based on /iSR. However, 
the field distribution is not consistent with the sharply 
defined stripe models that properly account for the elas- 
tic neutron scattering data and the precession of muons 
at faster time scales. 

To fit the Zeeman perturbed NQR within the stripe 
models, including the incommensurate spin density wave 
picture that is favored by /iSR, the observed lineshape re- 
quires a field distribution centered close to zero with very 
large distribution coupled with an increased quadrupole 
width. The magnitude of corresponding maximum Cu 
magnetic moments, ~ 0.15 /is, is suppressed relative to 
the '^0.3^B as determined by /LtSR studies. This, to- 
gether with somewhat narrower lineshape observed for 
Lai.g4Euo.2oSro.i6Cu04 and relatively high relaxation 
rate observed even at 350 mK, indicates that 

motional narrowing that is caused by remnant slow fluc- 
tuations which persist even at our base temperature. The 
presence of slow fluctuations has also been emphasized 
by Suh etal.M The fact that the 1.7 K and 350 mK spec- 
tra are of identical shape suggests that the decreasing 
frequency scale of charge and spin fluctuations has sat- 
urated and that there are many low lying excited states 
of the stripes, while the continuing NQR signal recovery 
may indicate that an increasing fraction of the sample is 
experiencing the saturated, low frequency fluctuations. 

The increased quadrupole width that is required to 
smooth the A and B lines into the low temperature 
lineshape indicate the presence of a broadly distributed 
charge environment which alters the electric held gra- 
dient site by site. Although we find in section [V that 
we can phenomenologically account for the ^'^Cu wipeout 
data based on slowing spin fluctuations, this quadrupole 
broadening indicates that spatial variation of the local 
charge environment does indeed exist throughout the 
sample. 



IV. SLOWING OF FLUCTUATION TIMESCALE 
FOR THE h PHASES 



11 



A. Preliminary considerations on the fluctuation 
frequency scales of stripes 

In general, the growth of short range antiferromagnetic 
spin-spin correlations causes the spin fluctuation fre- 
quency scale r — \/Tspin to gradually slow down. In the 
normal metallic state above in the high cuprates, 
the characteristic fluctuation frequency T of spin fluctua- 
tions decreases roughly in proportion to the temperature 
T. This peculiar slowing down is frequently referred to 
as "lj/T scaling," because the frequency dependence of 
x"(q, uj) satisfies a simple scaling law when the frequency 
Lo is namalized by the temperature T (or equivalently, 
by r).c3o One can define F as the peak frequency of 
the oj dependence of x"(q, w), or, more precisely, by fit- 
ting x"(q, w) to a relaxational function with appropriate 
frequencjjindependence. We fit x"(q, w) observed by Aep- 
pli etalxB for Lai.86Sro.i4Cu04 to a Lorentz oscillator 
form to deduce hT k, 8.8 meV at 35 K. By scaling the 
magnitude of hV obtained at 35 K for Lai.86Sro.i4Cu04 
linearly with temperature, we estimate hT « 17.6 meV 
at Tcharge = 70K. Nccdlcss to say, hT should depend on 
the hole concentration x (for exasnple, hT « 4meV at 
similar temperatures for x = 0.04E5). We expect a some- 
what smaller value of hT for x = 0.12 than for a; = 0.14 
because the dynamic spin-spin correlation length is some- 
what longer for x = 0.12 than for x = 0.15JZj However, 
since depends little on x for 0.12 < x < 0.20, we 

do not expect a major change in the magnitude of hT. 
We note that in Lai.8_2;Euo.2Sr2,Cu04 is nearly 

identical with that of La2-2:Sr2,Cu04. Hence, it is safe 
to assume that the spin dynamics above Tcharge depend 
little on Eu or Nd rare earth codoping. 

The important point to notice is that this mild T- 
linear slowing down of spin fluctuations with tempera- 
ture breaks down in the charge ordered segments of the 
Cu02 planes below TchargrT" 70 K. EPR measurements 
by Kataev and co-workersE£l indicate that spin fluctua- 
tions are slowing more quickly below Tcharge- Somewhat 
below the onset of charge ordering, spin fluctuations in 
some segments of the Lai.48Ndo.4Sro.i2Cu04 sample be- 
come slow enough at T^^^^"" = 50K to be considered 
static, and neutron scattering measurements begin to de- 
tect quasi-elastic scattering with excitation energy hu « 
1 meV. We emphasize that this does not mean that all 
the Cu spins have slowed down to hT f^i 1 meV. The fact 
that some Cu NQR signal is still observable even below 
Tcharge, and the ^^1/Ti measured for the obseryahlepSeg^. 
ments does not show any anomalous behaviorElEiJO'EEl 
indicates that the observable Cu NQR signals originate 
from other segments of the Cu02 plane where T is still 
smoothly slowing following the u/T scaling. One may 
imagine that phase separation sets in with microscopic 
length scales below Tcharge, and T begins to exhibit glassy 
slowing with a broad distribution because T slows down 
exponentially in the charge ordered segments. This is 
in sharp contrast with the narrow and well defined spin 



fluctuation frequency T in the ordinary metallic state of 
Cu02 planes above Tcharge- We also note that if charge 
order is a second order phase transition, the critical slow- 
ing of charge dynamics should exist even above Tcharge- 
In fact, ±iifi_resistivity data does show an upturn prior to 
Tcharge -Ej'O This mcans that precursors of glassy slow- 
ing may exist at temperatures slightly above Tcharge- 
Our earlier observation that the Gaussian component of 
^^1/T2 shows gradual crossover to Lorentzian due to mo- 
tional narrowing ~ 10-30 K above Tcharge supports this 
idea. 

The inelastic magnetic neutron scattering intensity 
reaches about 80% of the maximum value by 25 K, mean- 
ing that 80% of the spins are slowed to below the time 
scale of Tneutron ~ 10~^^ scc by 25 K. To our knowledge, 
no detailed inelastic neutron scattering data is available 
for La2-2,-y(Nd,Eu)j,Sra;Cu04 with x w ^ near Tcharge 
and Tspin other than a limited data set reported by Tran- 
quada etal.B Accordingly, the details of how the slowing 
of spin fluctuations deviate from the high temperature 
"uj/T scaling" towards hT « ImeV at 50 K as a func- 
tion of temperature is unknown. In what follows we 
smoothly extrapolate the low temperature behavior to 
the uj/T curve. 

In this temperature range, fiSTi measurements do not 
detect any static hyperfine fields. This means that Cu 
spins fluctuate so quickly that the hyperfine field is av- 
eraged to zero. With decreasing temperature, /xSR mea- 
surements begin to detect static hyperfine magnetic field 
that last longer than 0.1 /isec between 25 and 30 K. This 
means that the spin fluctuation frequency T has slowed 
down to '-^lO^ Hz in some segments of the Cu02 planes. 
At 30 K, the Cu NQR signal is almost completely wiped 
out as shown in Fig. |[ This indicates that normal metal- 
lic segments of Cu02 planes no longer exist below 30 K. 
By extending the oj/T scaling down to 30 K, we estimate 
the upper limit of the distribution of spin fluctuation fre- 
quency F ^ 10^^ Hz at 30 K. On the other hand, the 
lower limit from the slowest components of spins is ~10^ 
Hz. 

The line broadening of Zeeman perturbed NQR caused 
by static hyperfine magnetic fields takes place for 
^^^La NQR and ^^'^^Cu at an even lower temperature, 
T'jvqT"" « 18K and 8K, respectively (see Figs. | and 
^). This indicates that the fluctuation frequency spec- 
trum has slowed to a timescale comparable to the du- 
ration of a spin echo experiment. This time scale is set 
by the separation of two RF pulses, usually r = 12 /isec, 
i.e. if spins fluctuate much slower than hT w 10^ Hz, 
the hyperfine field would look entirely static and we ex- 
pect full hyperfine broadening. As explained in section 
III , our Zeeman perturbed NQR results seem to indicate 
that this condition is not quite satisfied in the present 
case even at 350 mK. In Fig. |l^ we present the summary 
of experimental information regarding the distribution of 
slowing spin fluctuations. 
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indirect spin contribution to wipe out can be triggerred by 
charge ordering. In section [V C , we analyze the sudden 
and dramatic increase of ^'^^l/TiT and its distribution 
that takes place below Tcharge based on the renormalized 
classical scaling, and deduce the spatial distr ibution of F 
as a function of temperature. In section IV D , wc analyze 
^^Cu and -'^^^La NQR wipeout effects together with mag- 
netic neutron scattering and /iSR data, again based on 
the renormalized classical scaling. The distribution of F 
deduced by this approach agrees well with that deduced 
from "^l/Ti. 



B. Spin stiffness in charge ordered segments and its 
relation to ^^Cu NQR wipeout 



Temperature [K] 

FIG. 10. The distributions of the spin fluctuation fre- 
quency, r, as a function of temperature for Nd and Eu 
codoped La2-iSr2,Cu04 for several values of x. For each mea- 
surement, the triangles, the small vertical line segments and 
the open circles correspond to the onset, the 50% point and 
the 80% point of the anomaly, respectively. In some cases data 
is not available and we estimate values from similar samples, 
as noted on the figure. The vertical positioning of the points 
indicates the approp riate timescale for each measurement, as 
discussed in section [VD, The thick dashed line across the 



The Cu nuclear spin-lattice relaxation rate "^l/Ti 
diverges exponentially in the-cenormalized classical scal- 
ing regime of both isotropicE3 and anisotropicO quasi- 
2d Heisenberg systems. As shown in Fig. |ll|, the expo- 
nential divergence of was successfully observed in 
undoped La2Cu04, and the fit to the renormalized clas- 
sical expjiessions allowed Imai etal. to estimate 27rps « 
1730K.lalia Since 2-Kps = 1.13J, this implies J « 1500K, 
in good agreement with estimatioH-based on neutron and 
Raman scattering measurements.ESEj 

Naturally, one would expect that measurements of 
^^1/Ti may allow us to determine thajSpin stiffness also 
in the striped phase below Tcharge -trEj Unfortunately, 
measured for the observable ^"^Cu NQR signal 
represents the spin fluctuation properties of the not yet 
stripe ordered segments (i.e. uj/T scaling), and exhibits 
littLs signature of exponential divergence even below 50 
KliilllZl where spin freezing is observed by neutron scat- 
tering. The absence of divergence in the measured value 
of can be understood based on the following con- 

siderations. First, we recall that, quite generally, nu- 
clear relaxation rates are inversely proportional to tha 
energy scales of spin fluctuations, such as 2'Kps and 
Since ^^1/Ti in hole doped high Tc cuprates, including 
the striped phase, exhibit nearly identical—values with 
undoped La2Cu04 at higher temperatures ,l3 the funda- 
mental energy scale of is set by the same bare 
J at high temperatures, even in the striped phase. On 
the other hand, from the measurements of the spin-spin 
correlation length, Tranquada et al. showed that the ef- 
fective spin stiffness 2Trp'^^f is as small as ~200K below 
Tcharge-" Wc notc that the fairly small spin stiffness can 
be considered a natural consequence of slowed charge dy- 
namics, because a hole intervening the Cu-Cu exchange 
path would locally diminish the value of exchange inter- 
In the following sections, we will demonstrate that the action J, and the effective spin stiffness 2'Kp''/ ^ is essen- 



top of each panel represents uj/T scaling where there is little 
distribution of F. The experimental data serves as anchoring 
points for the distribution of F that is used in the wipeout 
calculations, which is depictsdias a contour plot utilizing the 
solid black and gray curves.CJ The region between any two 
adjacent curves represents 12.5% of the volume of the sam- 
ple, with an additional 12.5% above the top curve and below 
the lowest. Dashed lines near F « 6 x 10^ Hz for x — 0.12 
represent the distribution F deduced in the Gaussian model 
as discussed in the text. The insets show the ^''Cu wipeout 
data and the simulated wipeout based on Eq. |l^. 



glassy slowing of spin fluctuations below Tcharge for x 



1 

8 

can be explained consistently based on the renormalized 
classical scaling of the Non-Linear-cr model w ith a re- 
duced value of spin stiffness 2-np\^^ . In section IV B, we 
estimate the numerical value 2-n-pl^-^ ~ 200 K based on a 
simple physical argument. We will also explain why the 



tially a spatial average of J at longer length scales. 

We now recall that all of the ^^Cu NQR signal in high 
Tc cuprates is near the detection limit of pulsed NQR ex- 
periments, and only a factor ~3 enhancement in 
and an accompanying change in ^^1/T2 is sufficient to 
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make signal detection impossible. That is why ^^Cu 
NQR measurepients in La2Cu04 could not be conducted 
below ~400 K.Ea The drastically reduced value of the spin 
stiffness in the charge ordered segments mean that once 
the charge order sets in, and ^^1/T2 in those seg- 

ments begin to blow up inversely proportional to I'Ep'i^ ^ . 
The fact that we do not observe a divergence of ^'^1/Ti for 
the observable parts of the NQR signal near Tcharge indi- 
cates that this divergence of ^^1/Ti in the charge ordered 
segments takes place very quickly. In fact, as summarized 
m Fig. 0, the Cu spin fluctuation frequency decreases 
by ~ 5-7 orders of magnitude between 70 K and 350 mK, 
roughly an order of magnitude every 10 K. This provides 
a natural explanation of why the nuclear resonance sig- 
nal from those segments experience dramatic and sudden 
wipeout due to the extremely fast relaxation times. This 
is the indire C | t s pi n contribution to wipeout triggered by 
charge order 
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FIG. 11. data as a function of tspiperature 

for undoped La2Cu04Mn (A), Lai.85Sro.i5Cu04Ej (x), and 
Lai.68Euo.2oSro.i2Cu04" (•). We also show ^^^l/Ti scaled by 
Eq. [n for the Eu codoped material (o). This quantity roughly 



til the low frequency spin fluctuations have slowed four 
more orders beyond the onset of ^^Cu wipeout. Accord- 
ingly, ^■^^La NQR is a better measure of spin 
fluctuations in the stripe ordered segments of the Cu02 
planes below Tcharge- However, we caution that ^"^^l/Ti 
is known to be dominated by processes other than Cu 
spin fluctuations above T ~ Tcharge, most likely electric 
quadrupole coupling with lattice vibrations, making in- 
terpretation of ^^^1/Ti above Tcharge difficult. 

In Fig. |ll| we plot both '^^l/Ti and ^^^1/Ti for 
Lai g8Euo.2Sro.i2Cu04. As initially discovered by Chou, 
Cho, Borsa and co-workers J2§E3 ^^^l/Ti shows roughly 
exponential divergence below T^vqa, in the underdoped 
regime of La2-a;Srj;Cu04. At a; = the ^^Cu data fits 
well to the renormalized classical form in the low T limit 
(T < 27rp,./2), written asB 
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0.35 (T/27rp,)3/2 



Zc Jha{l + T/2t:PsY 
he exp{2TTps/T) 



C = 0.27 



2Trpg 1 + T/2nps 



(8) 
(9) 



with the Q — (tt, tt) component of the hyperflne coupling 
constant Aq — — AB, the constant Zc = 1.18, and 
2ttPs = 1730 K. The spin wave velocity is c = \/2JaZc/h. 
We also found similar behavior of ^^^1/Ti in the striped 
phase of Lai.g_a;Euo.2Sr2;Cu04 and La2-a;Baa;Cu04 be- 
low Tcharge, whcrc the value of 2'Kp\^^ is greatly reduced. 
In fact, from Fig. 11 we can estimate 2-Kplff 



200 K, 

which is in good agreement with the previously meru 
tioncd result of neutron scattering, 2'Kp'^J^ — 200±50K.lI 
It is important to notice that the extrapolation of the low 
temperature renormalized classical fit to higher temper- 
atures, as shown by the dashed curves in Fig. |ll|, over- 
estimates 1/Ti observed for ^'^Cu NQR by an order of 



represents ^'-^ 1 /Ti of the unobservable ^^^Cu NQR signal in the magnitude above Tcha 



stripe ordered segments. The curves represent l/Ti in the 
renormalized classical scaling model with different values of 
as listed. Inset: the same data for the Eu codoped 
material on a linear scale for clarity. 



C. Spin fluctuation spectrum deduced from analysis 
of Ti recovery data 



Even though the wipeout of ^'^Cu NQR prevents us 
from observing the renormalized classical behavior in 
striped phase, one can utilize -'^^^La to probe the slowing 
of spin fluctuations. The hyperflne coupling between the 
Cu electron spins and ^^^La nuclear spins is two orders 
of magnitude smaller than for the ^^Cu site. Since 1/Ti 
is proportional to the square of the hyperflne coupling, 
this means that, if we ignore nonmagnetic contributions 
to the relaxation, ^'^^1/Ti at the La sites is four orders 
of magnitude less sensitive to the divergence caused by 
the slowing Cu spin dynamics as compared with ^'^1/Ti 
at the Cu site. 
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Thus 3JLa NQR signal does not suffer signal wipeout un- 



The renormalized classical expression for the nuclear 
spin-lattice relaxation rate 1/Ti (Eq. H) can be rewritten 
in terms of a frequency scale F = {c/S^)^jT /2'Kps as 



1 

— = 49 — 

Ti h'' F (1 



1 (r/27rp,)2 
- T/2Tips 



(10) 



which connects F to the observable quantity Ti, thus 
providing a natural route to obtain information on the 
fluctuation spectrum. Due to the strong wipeout of the 
^^Cu NQR signal, we are forced to rely on the "'^^^La re- 
laxation data, where wipeout is less severe and limited 
to the range 5 - 35K (see Fig. 13 for wipeout data). We 



(7) measure '^^1/Ti by the inversion-recovery method. The 
recovery of the nuclear magnetization may be flt to the 
solution to the standard rate equations. 
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for the I ± 



NQR transition where m(t) is 
the magnetization at a time t and 1/Ti is the single 
relaxation rate. This form is valid for magnetic relax- 
ation when the spectral function is constant across the 
three NQR transitions, which is generally the case for 
r ^ Jnqr- At 65 K, the fit of the recovery data with a 
single value of Ti is reasonably good, as seen in Fig. [T^a, 
but at 350 mK, the fit is quite poor. This led us to sim- 
ulate the recovery curves due to distributions of Ti . The 
distributions in 1 /Ti are displayed in the inset of Fig. [T^a, 
and the recovery is calculated by summing contributions 
from each segment of the sample, where the magneti- 
zation of each individual segment recovers per Eq. |l^. 
These simulated recovery curves are plotted as solid lines 
in Fig. [l^ a and provide a much better fit to the 350 mK 
data, while slightly improving the quality of the fit at 
65 K. The small distribution evident within the 65 K re- 
covery data does not necessarily reflect an intrinsic dis- 
tribution of the magnetic relaxation time Ti but may be 
due to extrinsic processes such as quadrupole relaxation. 
In general, the central value of the Ti distribution is very 
close to that which we find by force fitting to a single 
relaxation rate. As seen in Fig. [T^b, we found that the 
width w of the Ti distribution increases with decreasing 
temperature below Tcharge from w^20 to w~ 200-1000 
at 1.7 K and below. The local minimum near 15 K is 
caused by the fact that ^■^^La NQR signals with extremely 
fast ^^^1/Ti are effectively wiped out due to short T2 re- 
laxation rates (see Fig. |l^, hence the measured nuclear 
spin recovery does not include those contributions. In 
the absence of those contributions, the width w^lOO 
near 15 K should be considered as a lower bound on w. 

Using the RC form (Eq. |l^) with 27rp''/f = 200 K we 
can convert these distributions of 1/Ti to F without any 
additional parameters. For F 3> /nqr the conversion is 
straightforward; these points are shown in Fig. 12;. For 



F close to Jnqr the renormalized classical form may be 
rewritten as 
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{T/2nps 



r (i + /^Qfl/r2) (i + r/27rp,)2^ 



(12) 



where Jnqr = 18 MHz and we assume that the spin dy- 
namics follow the Lorentz form in the renormalized clas- 
sical regime. This function is now double valued in F, 
which, in addition to signal wipeout, makes the deter- 
mination of F in the region 5-35 K difficult. In that re- 
gion, F ~ Jnqr, which implies fast relaxation and signal 
wipeout. From Ti measurements alone it is impossible 
to know how much of the distribution of F is greater or 
less than /nqr- All we can say for certain is that experi- 
mentally the value of 1/Ti is underestimated because the 
fastest sections are not observable. However, at 350 mK 
these difficulties disappear because all of the signal in- 
tensity is recovered and we know that F < /nqr- Wc 



can now safely employ Eq. to determine the distri- 
bution of spin fluctuation frequencies, which is centered 
at 3 X 10^ Hz with a width of more than two orders of 
magnitude. 
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FIG. 12. 



'La NQR recovery data at 350 mK and 65 K 



for Lai.68Euo.2oSro.i2Cu04 fit to single values of Ti (dashed 
curves) apd to distributions of Ti (solid curves) as given in 
the inset.Cil The width w of the distributions (defined as the 
quotient of the two Ti values at the half maximum positions of 
the distributions) are plotted as a function of temperature in 
(b) for Lai.875Bao.i25Cu04 (A) and Lai.68Euo.2oSro.i2Cu04 
(•). We note that the width is underestimated during NQR 
signal wipeout (the shaded region) and with this in mind we 
draw small vertical arrows from the most affected measure- 
ments and a dotted curve as a guide for the eye. In (c) the 
fluctuation frequency F deduced from Ti (from forced fit with 
single Ti, (x); central value of best fit distribution (o) with 
F WHM as given by the vertical bars) . At 350 mK we display 
data from a Lorentzian extension (o) of the form of Ti and 
from a Gaussian extension (o) (central value and FWHM po- 
sitions marked with matching symbols for clarity). The solid 
and dashed curves are reproduced from Fig. ^o]b and act as 
guides for the eye and the dotted curve represents the form 
exp(-27rpf ^/T) with 27rpf = 200 K. 



In the preceeding analysis for F < Jnqr we chose to 
extend the RC form of l/Ti with a factor of 1/(1 -I- 
InqrJ^'^) (l^l- 12)- This assumption is certaintly jus- 
tified for the paramagpetic spin dynamics above 70 K 
where the neutron datao fits well to a Lorentzian. How- 
ever, there is no reason to believe that the very slow dy- 
namics in spin stripes at 350 mK have to be governed by 
the same type of relaxational mechanism. For instance. 
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if scattering by meandering hole rivers takes over as the 
dominant mechanism of spin fluctuations at low temper- 
atures, the form of the relaxation function may change. 
Let us try another sensible choico—the form based on 
Moriya's Gaussian approximation)^ which would give 
frequency dependence in 1/Ti as 
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(13) 



We also employ the form to deduce F at 350 mK, and 
find a much narrower distribution centered at a higher 
frequency, 5 x 10^ Hz. We will discuss the implications of 
choosing either the Lorentzian or Gaussian from of l/Ti 
in section IV E| . We note that in simulating the Ti recov- 
ery curve we continue to employ Eq. |ll| even for T close to 
and less than Jnqr- In that region the spectral function 
is not flat over the three resonant transitions, invalidating 
the basis of Eq. Continued use of Eq. |ll| is justifiable 
because the distribution in extends over several 

orders of magnitude and should be much greater than 
the potential error induced through the continued use of 
Eq. |ll|. 



D. Numerical simulations of NQR signal wipeout 
and recovery, neutron, and /iSR 

As previously defined, wipeout is the amomalous loss 
of NQR signal intensity. For La2_2;_3^(Nd,Eu)j;|Srj;Cu04 
with 0.12 w i, the onset temperature TjvQflEd of the 
^^Cu NQR wipeout agrees very well with the onset tem- 
perature. |TLi__,e of short-range charge order detecteiLhy. 
neutronHCnalHj and x-ray scattering measurements£II'E£l 
Moreover, the fraction of wipeout F{T) shows a charac- 
teristic temperature dependence which is nearly identi- 
cal to the charge order parameter observed by scattering 
techniques (i.e. the square root of the scattering intensity 
of the charge order peak), leading us to argue that the 
temperature dependence of the wipeout fraction F(T) 
is a good yneasuralj of the charge order parameter. Re- 
cent i^^La NMR studies of La2Ni04+5 by Abu-Shiekah 
et al£3 also showed that the temperature dependence of 
the wipeout fraction of ^^^La NMR in a high magnetic 
field (with an onset temperature of ^ 200 K) agrees quite 
well with the charge stripe order parameteiH observed 
by elastic Bei,itron scattering by Tranquada, Wochner and 
coworkers EjQ 

Naively thinking, the slowing charge fluctuations can 
make direct contributions to wipeout through extreme 
NQR line broadenina., similar to the case of conven- 
tional CDW systems, E3l3 and also, through extremely 
fast quadrupole relaxation times.OcS It is difficult to ac- 
cess how strongly the charge fluctuations near Tcharge 
affect the NQR properties of the missing signals because 
the signals from the already stripe ordered segments are 
not observable, similar to tho-all or nothing behavior of 
wipeout observed for Cu:Mn.E3 



In addition to the direct charge contribution to wipe- 
out, we have also pointed out that charge order turns on 
low frequency spin fluctuations below Tnqh within the 
charge ordered domains, which in turn causes wipeout 
by making ^^Ti and '^^2 too short to detect the ^^Cu 
NQR spin echo signal.tj That is, even if charge fluctu- 
ations are not strong enough to cause the wipeout of 
^^Cu NQR directly, charge order creates patches of the 
Cu02 plane that behave similarly to conventional spin 
glasses at slow time scales. As already noted by Huat 
et. ai, dilute magnetic moments in metals (e.g. Cu:FeL2l) 
are known to cause wipeout of the NMR signal of the 
host metal. Restating this, we can say that the slow- 
ing charge dynamics are the indirect cause of ^'^Cu NQR 
wipeout by triggering the drastic slowing of spin fluctu- 
ations within the regions where holes are trying to form 
charge rivers. Unfortunately, these important statements 
from Hunt et. al. and Singer et. al regarding the effects of 
strong spin dynamics on wipeout have been overlooked 
by some readers, including the authors of a recent publi- 
cation from Los Alamos.Ej 

Hunt et al. cited several pieces of experimental evi- 
dence that are consistent with the picture that anomalous 
slowing of spin dynamics are caused by slowing charge 
dynamics. First, inelastic neutron scattering measure- 
ments conducted at low energy transfers (1.75-3.5 meV) 
seem to show an anomalous enhancement of the low fre- 
quency spin fluctuations near and below Tcharge, even be- 
fore quasi-elastic neutron scattering with energy transfer 
less than 1 meV detects frozen spins.tl 

Second, the spin-lattice relaxation rate ^^^1/Ti at the 
^^^La sites begins to exhibit an anomalous upturn be- 
low Tmqr. This enhancement of ^'^^1/Ti is followed 
by a roughly exponential divergence at lower temper- 
atures, i39i/Ti ~ exp{2T: p"/ f /kBT), where 2np''Jf is 
the effective spin stiffness which is of the order 10^- 
lO^K, as already shown by earlier comprehensive ^■^^La 
NQR measurements by Chou, Cho, Borsa. Johnston and 
co-workers in underdoped La2_a:Sra;Cu04E3'EZl Recalling 
that the exponential divergence of 1 /Ti is a typical signa- 
ture of the low temperature renormalized classical scal- 
ing behavior oLB.. — ^ quantum Heisenberg svstems in 
both isotropicOEj and anisotropic 2d systemsjZ3 the ex- 
ponential divergence of below T^qj^ is consis- 
tent with a rather drastic reduction of the effective spin 
stiffness /Og-^-^ caused by slowing charge dynamics. The 
spin-spin correlation length measured by Tranquada et al. 
also supports the, renormalized, small spin stiffness with 
2np'g^^ « 200 K.H In contrast, we recall that the funda- 
mental energy scale in hole doped Laa— ™Sra:Cu04 is set 
by 2t: ps Rj— 1730K as shown by NQROO and neutron 
scattering.E3 High energy scales of the order of 600 K are 
also inferred from a scaling_analysis of the unifopn spin 
susceptibility by JohnstoiO and Hwang etal..L3 This 
suggests that drastic reduction in spin fluctuation energy 
scales from the order of 10^ K to 10^ K occurs near the 
charge ordering temperature. This has an important con- 



16 



sequence on the wipeout phenomenon, because quite gen- 
erally the overall magnitude of nuclear relaxation rates 
and ^^1/T2 are inversely proportional jta the spin 
fluctuation energy scale ^H/Ti,2 oc I / 2t: p^J f H The ob- 
served relaxation rate ^^\/Ti ^ 2000 sec^^ is quite fast 
even above Tcharge in hole-doped high Tc cuprates, while 
^^Cu NQR signals become unobservable when relaxation 
rates exceed ~ 6000 see" ^ for measuremeiits with 

the quantization axis along the crystalline c-axis.E3 When 
exceeds eOOOsec"^ the relaxation times ^^Ti 
and ^^T2 become shorter than the limit of experimentally 
accessible times. In other words, the ^'^Cu NQR signal 
in all high Tc cuprates is on the verge of wipeout due to 
strong low frequency spin fluctuations. This means that 
a factor 2~3 reduction in spin stiffness p^-^-^ near Tcharge 
caused by slowed charge dynamics would immediately 
wipe out ^■^Cu NQR signals. This explains why the ^^Cu 
NQR wipeout fraction accurately tracks the charge order 
parameter as measured by scattering experiments if one 
interprets the latter as a measure of the volume fraction 
of charge stripe ordered segments. 

Third, the Gaussian component of spin-spin relaxation 
process ^^•^^1/T2 at the ^^'^^Cu sites crosses over to 
Lorentzian (single-exponential) decay below Tcharge by 
motional narrowing, effects due to fluctuating magnetic 
hyperfine fields .E2rEJ The dominance of spin rather than 
charge fluctuations to the mechanism of the Lorentzian 
spin echo deca5£||«a,s suggested by the isotope ratio be- 
tween ^'^'^^Cu.IHJ'EJ This crossover of spin echo decay, ac- 
companied by NQR/NMR intensity wipeout, is known 
to be a typical signature of slowing dassp. of conven- 
tional spin-glass systems such as Cu:Mn.E3o Thus all of 
these rcsuJls as well as the dramatic increase in the EPR 
linewidthllj are consistent with a simple picture that spin- 
glass like behavior sets in when charge dynamics slow 
down, and the combined effects of slowing charge and 
spin stripe fluctuations can naturally provide a qualita- 
tive account of the wipeout behavior for x > ^ below 

TnQR ~ Tcharge- i—i 

In a recent article, Curro, Hammel and coworkerala 
made an important observation that enabled quanti- 
tative estiiiiation of the indirect spin contribution to 
wipeout JIJ'EJ The essence of the analysis of Curro and 
Hammel is as follows. First, instead of saying-tliat ex- 
tremely short ^^Ti and ^^T2 result in mpeoutEB'tiJ they 
note that the ^^Cu NQR signal cannot be detected if^^Ti 
and ^^T2 are shorter than certain cutoff values labeled 
Q^rpcriticai 63 ycritica/ ^ Using Several assumptions, 
they estimate these cutoff values. Qace 632^'"'*"=°' is de- 
duced, Curro etal. use the BPP formcjof Ti to determine 
the corresponding spin fluctuation frequency 63pi«zpeout^ 
The wipeout fraction is the weight of the spin fluctua- 
tion distribution that is slower than 63pt«ipeoMt bg^ause 
the slower segments have ^^Ti and ^^T2 that are too fast 
to allow the detection of the spin echo. 

Inspired by their analysis, we attempt to fit not only 
the wipeout of the '^^Cu NQR signal below Tnqr « 70K 
but also its recovery below 8 K, the wipeout of ^'^^La NQR 



signal below 40 K and its recovery below 18 K, the spin 
order parameter measured by elastic neutron scattering 
below 50 K and the pSK asymetry below 30 K by tak- 
ing into account the qualitative change in spin dynamics 

below Tcharge- 

We take the following approach, which relies pri- 
marily on experimental data instead of theoretical as- 
sumptions. In order to set 63pcMto// note that 
s^l/Ti w 2000 sec" 1 at 70 K, at which temperature 
we know that the energy scale of spin fluctuations is 
hV w 15 meV from neutron scattering. Fro«i earlier 
^■^Cu NQR studies in paramagnetic La2Cu04j^ we em- 
pirically know that the loss of ^^Cu NQR intensity due to 
fast relaxation occurs when the relaxation rate exceeds 
e^l/T^riUcai _ 6000sec-i. Following Eq. |l|, this implies 
that 63ptDipeout approximately a factor of three lower 
than the value of hV w 15meV at 70 K. In this way, we 
deflne esp^ipeo^t ^ 5 j^^^y^ ^YyicIy is 1.2 x lO^^ Hz (so that 
l/63p™<peo«t _ 0.83 psec). 

The cutoff for the wipeout of ^'^^La NQR signal is found 
in an analogous manner. From our experments shown be- 
low in Fig. 16 and Fig. |l^, we see that ^^^1/Ti rises to ap- 
proximately 80 sec" ^ when wipeout is most severe. The 
remaining portion of the sample is on the verge of wip- 
ing out and has Ti very close to i39 j^/ycntica/^ Thus we 
take 139 j^/j'crjtica/ ^ gOscc"^. We now need to convert 
this value of ^'^^ 1 /T^^ztzcai ^j^g corresponding frequency 
scale. Using the renormalized classical form (Eq. in 
the high frequency limit (T ^ ujn) we see that 
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(14) 



where wipeout occurs at 70 K for ^^Cu and at 40 K for 
^^^La, and we take 27rps w 200 K as determined above. 
Since we know the ratio of the gyromagnetic ratios to be 
In/^^ln = 0.533 and the ratio of the hyperfine fields at 
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the two nuclear positions to be '^■^Hhf/ Hhf = 0.012, as. 
inferred from Zeeman perturbed NQR in La2Cu04j23'L£l 
we immediately calculate i39r^"'P'=™* = 5.1 x 10^ Hz. We 
note that the ^^^F{T) data is taken with a fixed delay 
time T = 40 /isec, which is justified due to the slow relax- 
ation rate ^•^^1/T2 at the La site (see the inset of Fig. p^ . 

At low temperatures, the ^'^Cu and ^'^^La NQR signal 
that was wiped out, begins to reemerge. The return of 
signal begins at about 8K for ^'^Cu and 18 K for ^^^La. 
Within our model, the reemergence of signal is not sur- 
prising, in fact, it is to be expected if the frequency range 
falls low enough so that Ti and T2 become longer. Us- 
ing the Lorentzian extension of the renormalized classi- 
cal form (Eq. |l^), with the NQR resonance frequency 
fNQR taken as 6 MHz and 35 MHz for ^^^La. and ^^Gn, 
respectively, that the relaxation rate 1/Ti achieves its 
maximum value when the characteristic fluctuation fre- 
quency r coincides with fNQR- Below that, additional 
slowing of the fluctuations (decreasing of F) reduces the 
relaxation rate 1 /Ti , so that at some low frequency scale, 
which is defined as r»'ecot)erj/^ ^^le following condition is 
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satisfied: i/T^{T'-<"""'<''-y) = i/T-^{Y^^peout^ ^ ^jr^iUcai 
where we include temperature factors from Eq. We 
note that r'''^™'"^''?' < /„ < r'"*'"^™*. Within this model, 
it is straightforward to deduce that that i39r'~'^™''ery ^ 
2.9 X 10^ Hz and earr-ecoi-ery ^ 4 X 10* Hz. This means 
that the slowest components within that sample slow to 
2.9 X lO'^ Hz at approximately 18 K, providing ^^SLa NQR 
signal recovery, and continue to slow through 4.0 x 10* Hz 
at 8K where they are manifested as recovered ^^Cu NQR 
signal. 



possible contribution of the slowing of charge dynamics 
to Cu NQR wipeout effects. 

The simulated wipeout fraction F{T) is the fraction of 
the NQR signal that is not observable, which corresponds 
to the portion which experiences fluctuations that are 
slower than T"'*?"^™* but faster than r''^™^''''?^ . This can 
be written as, 
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FIG. 13. Four sets of experimental data and correspond- 
ing numerical simulations based on the model described in 
the text, using the L»(r, /) shown in Fig. (a) "^^Cu wipe- 
out in Lai.68Euo.2Sro.i2Cu04, (b) the ^"^^La wipeout from the 
I ± I) <-» I ± I) transition in Lai.68Euo.2Sro.i2Cu04 (The 
wipeout behavior at s « i is nearly identical for Eu, Nd, 
and Ba co-doped samples.) (c) the spin order parameter in 
Lai.48Ndo.4Sro.i2Cu04 deduced from elastic neutron scatter- 
ing as the sftuare root of the scattered intensity, from Tran- 
quada et aLp and (d) the /xSR asymetry factor A_Y-/yAtotai x 
in Lai.475Ndo.4Sro.i25Cu04, from Nachumi etal.E3 

Combining all of these pieces of information regarding 
the spin fluctuation time scales in the stripe phase, we 
construct a chart showing the anchoring points of the 
temperature dependence of spin fluctuation frequency 
scale r, as presented in Fig. [l^b. Using these anchor- 
ing points, we smoothly interpolate to define a single 
distribution function D(T, /), which gives the relative 
fraction of the sample volume that experiences fluctua- 
tions at frequency / at temperature T. This distribution 
is depicted by the black and gray contours in Fig. |l^. We 
begin with an educated guess at the form of distribution 
function, and reflne the function as we iteratively pro- 
duce a better and better fit to the experimental data as 
presented throughout this section. In this way we deduce 
the actual fluctuation spectrum based on the experimen- 
tal data. Since there are no experimental constraints for 
r in the range above 63p»*peo«t jggg ^^^^ ^jj, 

scaling line, we are forced to extrapolate in that region. 
In Fig. |l^ we depict smooth extrapolations of D{T, /) 
that terminate at the lo/T scaling line at a temperature 
above Tcharge- Howcvcr, we note that we are ignoring the 



where the denominator acts as a normalization condition 
on D{T, /) at every temperature. For the neutron scat- 
tering and /iSR measurements, the fraction of the sample 
that has ordered is given by a similar form. 



Jl^"^"^ D{TJ)df 
S^D{TJ)df 



(16) 



where probe is either ^SR or neutron. 

All four of these quantities {^^FiT), '^^^F{T), 
G"^"*™"(T) and G''^'^(r)) have been numerically com- 
puted over a range of temperatures using the distribution 
D{T, f) presented in Fig. |l^, and the resulting curves 
are plotted in Fig. [l^ along with the relevant experimen- 
tal data. The simulations work remarkably well in repro- 
ducing both the ^^Cu and ^^^La wipeout data, as well as 
the neutron and /xSR results. As mentioned previously, 
the quality of these fits is the result of iterative alterations 
of trial distribution functions that eventually led us to the 
distribution D{T, /) shown in Fig. p^. We find that a 
distribution function based on a Gaussian in log/ with 
different widths to higher and lower frequency works well. 
However, the use of an iterative method means that we 
cannot exclude the possibility that a somewhat different 
/ and T dependence of D{T, /) may work as well as this 
one, although it is difficult to imagine that a qualitatively 
different form would succeed, given the large number of 
experimental constraints that must be satisfied. 

Although we do not explicitly utilize the nuclear spin 
lattice relaxation rate 1/Ti in the analysis of this section, 
our choice of the form of Ti for F below the Larmor fre- 
quency does play an important role in setting the vertical 
position of the anchoring points in Fig. |lO|. T o this point 

^qll). 



we have utilized t he Lo rcntzian form (Eq. |12|), but as dis- 



cussed in section IV C , we might do better to employ the 
Gaussian form (Eq. |13|). Although this change does not 
alter any of the analysis for F > 5 x 10^ Hz, it does sub- 
stantially change the behavior for smaller F. Specifically, 
the cutoff frequencies for the recovery of both ^'^Cu and 
^^^La NQR signal are shifted upwards to ^ 5 x 10^ Hz, 
causing the distribution D{T, f) to saturate at low tem- 
peratures as shown by the dashed curves in Fig. ^oJd and 
in Fig 12. We will discuss the merits of the Gaussian 



form of 1/Ti in section IV E 
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FIG. 14. Slices of the frequency distr Aution D{T, f) as de- 

I Lj The shaded rectan- 
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duced from the analysis of section 
gles represent the spectral weight in the sharply defined uj/T 
energy scale (70 K is the right most, 50 K in the center and 
30 K at the left), which decreases in magnitude at lower tem- 
perature as spectral weight is shifted to the contribution from 
the stripe ordered regions represented by the broad curves. 
The dashed curves arise from a Lorentzian extension of the 
form of Ti (Eq. |l^) and the shaded curve corresponds to a 
Gaussian extension (Eq. at 350 mK. 



E. Summary of fluctuation timescales and 
comparisons to related works 

Through an analysis of both the 1 /Ti recovery data 
and the ^^Cu and ^'^^La wipeout data we arrived at two 
estimates of the spin fluctuation spectra as a function 
of temperature as shown in Figs. 10) 3 an d As an 
alternative presentation of the data in Fig. |10^, we con- 
struct Fig. which contains slices through the distri- 
bution D{T, f) at selected temperatures. We find that 
there is a sharp change from the high temperature tu/T 
dominated behavior to the much faster exponential slow- 
ing of spin fluctuations below Tcharge- We understand 
this qualitative change in the magnetic fluctuation spec- 
trum to be due to charge ordering. The two methods 
of deducing F are in qualitative agreement, even though 
the quadrupole contribution to relaxation and wipeout 
makes the approach based on ^'^^1/Ti somewhat ques- 
tionable. The differences between the two methods are 
smaller than the effects of choosing a Lorentzian or Gaus- 
sian extension of the form of Ti. Without additional 
information we would be able to proceed no farther — 
however we know from the lineshape analysis that our 
350 mK lineshape is still experiencing some motional av- 
eraging. In the Lorentzian schemes more than 50% of 
the distribution of F is lower than F = 3 x 10''' Hz, which 
corresponds to SO^sec, the duration of a pulsed NQR 
experiment. This indicates that motional averaging is 
unlikely to be effective for the bulk of the sample and we 
would expect to see frozen structures in the low tempera- 
ture Zeeman perturbed spectra, which are not observed. 
However, using a Gaussian model we find that almost all 
of the sample experiences fluctuations faster than the ex- 
perimental duration, indicating the likelyhood of strong 



motional averaging. Thus, reality is closer to the Gaus- 
sian form in order to preserve motional averaging. We 
also note that the low temperature saturation of F that 
naturally arises from the aforementioned analysis within 
the Gaussian form is consistent with the saturatioix. of 
the spin-spin correlation length observed below 20K.EI 

Before continuing, we point out some qualitative simi- 
larity between D{T, /) for the present case of the striped 
cuprates below Tcharge and for the conventional spin 
glass system Cu:Mn investigated by Murani by cold neu- 
tron scattering studies.u It is also important to note 
that at the slow timescales of NQR, the striped phase 
never looks to be striped, but rather behaves more or 
less like a conventional spin glass. This similarity with 
the cluster spin glass phase h aa bge ji repeatedly empha- 
sized by Borsa and coworkersrnrTLj since before the dis- 
covery that the spjipi-iglass phase is actually striped by 
neutron scattcringfl^El However, we note that the striped 
phase does differ from conventional spin glasses, even 
at slow timescales. Specifically, we point out that the 
real part of the susceptibility shows only small Curie 
enhancement, indicating the scarcity of free spins. In 
fact, comparing Lai.87Bao.i3CuCl£3 and Zn impurity 
doped Lai.85Sro.i5Cuo.98Zno.o204ca where 2% of the mo- 
ments are intentionally made to be free, we find that 
the Curie component of the susceptibility is nearly two 
orders of magnitude smaller in the Ba doped material 
below Tcharge than in the Zn doped case. This implies 
that the number of free spins in Lai.87Bao.i3Cu04 is very 
small. It is more likely that the spin glass like behavior is 
caused by slow fluctuations of self organized, short range 
ordered segments with spin and charge density modula- 
tion. 

Recently, Curro, Hammel etaLo and Teitel'baum 
et alE3 have reported methods of simulating the wipe- 
out effects. In both cases the authors attribute ^^Cu 
NQR signal wipeout to slowing spin fluctuations. We 
emphasize that the importance of slowing spin fluctu- 
ations on NQR signal wipeout through extremely large 
relaxation rates ^'^1/Ti and ^^1/T2 (the mdirect contribu- 
tion), has already been discussed repeatedly in our ear- 
lier publications .E2Hlil The contribution of Curro, Ham- 
mel and coworkers is in introducing of a numerical cut- 
off to fit_wipeout effects, making quantitative analysis 
possible.Ej Their pioneering work opened a new path to 
gain quantitative information regarding the distribution 
of the spin fluctuation frequency F. This is an important 
idea that is applicable to the analysis of various NMR 
data and they should be fully credited for their contribu- 
tion. However, Curro et al. overlooked the fact that the 
importance of spin fluctuations to wipeout (the indirect 
contribution) had originally been proposed by us, and 
quoted our publications only in the context of slowing 
charge dynamics. We believe that this inaccurate quota- 
tion is misleading at best. Even though they have implied 
that their way of understanding the wipeout effects was 
entirely novel, we must emphasize that the physical pro- 
cesses utilized in their analysis have little difference from 
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the ideas outhned in Hunt et al. regarding the indirect 
process, except for one crucial point. We believe that 
the onset of charge order (or slowing of charge-dynam- 
ics for 2: < i due to random localization effectsU) is the 
cause of the drastic changes in the magnetic properties 
that result in NQR signal wipeout. On the other hand, 
in the analysis of Curro and coworkers it is assumed that 
glassy slowing of Cu spin fluctuations continuously pro- 
ceeds from 300 K to 4 K without any qualitative changes 
in spin dynamics through Tcharge- That is, Curro etal. 
assume that the same activation temperature dependence 
of r(T) c>c exp{—Ea/kBT) properly represents the slow- 
ing of Cu spin fluctuations both above and below Tcharge- 
They fit below 30 K (< Tcharge) to deduce the 

activation energy Ea and its distribution A. Ignoring 
the drastic change in the charge degree of freedom near 
Tcharge and the resulting changes in spin dynamics as 
reflected in EPRJIS etc., these parameters are 

taken to remain unchanged even above Tcharge- 

We note several inconsistencies arising from the as- 
sumptions of the Los Alamos group. We point out that 
the fit of ^^^l/Ti by Curro and coworkers fails on their 
own data set above 30 K. In addition, within the model of 
Curro et ai, ^^1/Ti has the same Gaussian distribution as 
r(r) with the same activation type temperature depen- 
dence, c>c exp{Ea/kBT) which is taken to hold for 
all temperatures. Unfortunately, this assumed exponen- 
tial increase of ^^1/Ti is in sharp contradiction with the 
experimental finding that decreases with decreas- 
ing temperature down to Tcharge as shown in Fig. We 
note that the assumed exponential form above Tcharge is 
essentially equivalent to extrapolating along the 
dotted lines of Fig. 11 between Tcharge and 300K. Thus, 



there is no justification for the assumed forms of r(r) 
and above Tcharge- Since the calculated wipeout 

fraction by Curro et al. relies on the exponential temper- 
ature dependence of ^'^l/Ti, we must conclude that qual- 
itative success in reproducing wipeout starting at as high 
as 150 K is merely a coincidence that depends strongly 
on their choice of 63p«.»peout^ 

Finally, although we believe that slowing of charge dy- 
namics is of extreme importance and that the analysis 
of the Los Alamos group is not apphcable for a: > |, we 
note that their assumed form of r(T) is similar to our 
finding for x = 0.07 below TpfQu. However, we note that 
in this region the beginning of charge localisation has 



been obs erved by Ichikawa et 
tion 



I at - T 



NQRi 



(See sec- 



V C for a discussion of the underdoped materials. 



Hence, one may consider the slowing of charge dynamics 
as the cause of the glassy slowing of spin fiuctuations re- 
gardless of spiatial coherence, which is gradually lost with 
decreasing x.tHlEd 

After atL, early version of this work was presented 
elsewhere,E3 the NQR wipeoujUjias also been simulated by 
a group at Leiden UniversityCJ who have independently 
pointed out that ^'^^La NQR intensity provides valuable 
information regarding F. They base their analysis on 
^^^La spin-lattice relaxation data in Nd codoped materi- 



als. Although we agree that ^^^1/Ti is a good place to 
initiate a study of the fluctuation spectrum and Cu wipe- 
out, we note that nuclear relaxation in the Nd codoped 
materials is usually|-dpminated by the large jMd moment, 
as in NdBasCuaOyN and Nd2_^Ce^Cu04.N This also 
seems to be the case in Nd doped La2_a;Sra;Cu04 as 
shown in Fig. |l^. Furthermore, we find that ^^^1/Ti 
in the Nd codoped materials is roughly independent of 
hole doping.ll23 Bearing this in mind, it seems that their 
analysis would then predict identical wipeout for differ- 
ent hole concentrations x, which is certaintly not the case 
(see Fig. IH). 



V. Nd AND Eu CODOPED MATERIALS AWAY 
FROM X PS 1 

A. Construction of a unified phase diagram 

Now that we have reached a good quantitative un- 
derstanding of the spin stripe fiuctuations and its tem- 
perature dependence for the magic hole concentration 
of a; « ^, we seek a comprehensive understanding of 
stripe fluctuations across the entire phase diagram of 
Lai.6-2;Ndo.4Sra;Cu04 and Lai.8~a;Euo.2Sr2;Cu04. The 
temperature dependence of the ^'^Cu NQR wipeout 
^^F{T), i39La NQR intensity, and ^^^1/TiT are pre- 
sented in Figs. |l^, |l^, and |l^, respectively. Utilizing 
this information, as well as available elastic neutron scat- 
tering and /i SR measurements, we follow the analysis of 
secton [VD| to deduce the temperature dependence of 
the spin fluctuation spectrum F for hole concentrations 
of X = 0.07, 0.12, 0.16 and 0.20. The results are depicted 
in Fig. |l^ together with the simulated ^'^Cu wipeout. 

Another way to illustrate the temperature dependence 
of the spin dynamics of stripes is to plot contours of char- 
acteristic frequency scales on the temperature (T) versus 
hole concentration (x) phase diagram. In Fig. |l|, we 
collect data from various experiments and plot the tem- 
perature at which each experimental probe experiences 
an anomaly. The anomaly occurs at the temperature at 
which the spin fluctuations of some parts of the CUO2 
plane slow to the inherent frequency of that probe. Us- 
ing the anomaly temperatures as a boundary, we section 
off areas of the phase diagram into various shades of grey. 

The white area at the top of the plot corresponds to 
the temperature region with very fast spin fluctuations, 
where stripes can be considered to be completely dy- 
namic. In other words, this is the region where the so- 
called Lu/T scaling holds for spin dynamics, and there is 
a well-defined spin fiuctuation time scale whose distribu- 
tion is minimal. (It should be noted that there is NMR 
evidence for some spatial inhomoginiety even in a similar 
high temperature regime, as recently reported by Haase, 
Slichter, and co-workers on the basis of some ingenious 
NMR measurements and analysis,t22l and by a report-of 
the frequency dependence of ^^1/Ti by Fujiyama etal. .12^ 
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FIG. 15. (a) The fraction of wiped out signal, '^^F{T), in 
Lai.6_i:Ndo.4Sra;Cu04. The low temperature structural pha^ 
transition temperature Tltt at x « | is Tltt ~ 70K£3 
The lines are guides for the eye and only data points 
below Tnqr are plotted. (b) is the same pl^t but for 
Lai.8-i:Euo.2Sri:Cu04. At a; « |, 
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16. Plot of the 
Euo.2Sri,Cu04 for the 



'La NQR signal intensity in 
±|) <-> |±|) transition at a fre- 
quency 3uQ « 18MHz {x = 0.07, o; 0.12, A; 0.16, V; 0.20, x). 
The intensity was measured for a fixed pulse separation time 
of r = 40 /xsec. In the inset we show the spin echo decay for 
Lai.6sEuo.2oSro.i2Cu04 at different temperatures. All curves 
are guides for the eye. 

The darkest grey region shows where the spin fluctua- 
tions are so slow that the hyperfine field is nearly static 
even at the NQR times scales. This region is obtained 
based on the appearance of hyperfine broadening in the 
recovered NQR signals of ^'^Cu. Between the lightest 
and darkest grey areas, which correspond to the wipeout 
and recovery of ^^Cu NQR signal, there are five separate 
regions depicted in Fig. |l^. The boundary marked by 
solid circles is set by the temperature at which ^^^1/TiT 
crosses the value O.OSsec^^K^^ in Lai.8-a;Euo.2Srj:Cu04 



(see Fig. is somewhat larger at high tem- 

peratures for X = 0.20, so in this case we use the temper- 
ature (25 K) at which the enhancement of be- 
gins. The next boundary is set by Tspin, the onset of the 
elastic neutron scattering from spin stripes correspond- 
ing to a frequency scale of 0.5 — 1.0 meV, based on the 
energy integration windows of those experiments. The 
next area marks the onset temperature of ^'^^La NQR 
wipeout. Open diamonds mark the onset temperature 
for muon spin rotation T^si?. arising fmm static order on 
the timescales of /iSR measurements. t3 At temperatures 
lower than 20 K the NQR signal begins to recover for the 
^^^La nuclei. These temperatures are marked by open, 
downward pointing triangles in Fig. |l^. Also shown in 
Fig. ^ are the onset temperatures Tcharge for shortojange 
charge order as nhtained from neutron scatteringE3 and 
x-ray scattering JlZllij the transition to the superconduct- 
ing phase in Lai.6-a;Ndo.4Sr2;Cu04 (a white line), and 
a dotted black line (as a guide for the eye) that follows 
Tcharge and the Cu wipcout inflection points. The im- 
portance of the inflection points will be explained below. 




20 40 611 KO 

Temperature [K] 

FIG. 17. (a) Plot of ^'^^l/TiT of the | ± |) ^ | ± |) tran- 
sition for La site in Lai,g_a;Euo.2Sra;Cu04 where the symbols 
used for each hole concentration x is listed in the figure. The 
horizontal line marks the constant ^^"^l/TiT ^ 0.05 sec~^K~^ 
value and the arrows mark Tnqr for each material x = 0.20, 
0.16 and 0.12 from left to right, (b) The same quantity, for 
several samples with hole concentration a:: « | . Note that the 
Nd codoped system has an enhanced relaxation rate due to 
the Nd moments. The bump in the data for the Ba doped 
sample near 60 K arises from the LTT structural transition. 
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FIG. 18. Phase diagram for Lai.6-i:Ndo.4Sra;Cu04 and 
Lai.8-2!Euo.2Sri;Cu04, showing the upturn in ab-plane Be-i 
sistivity temperature Tu found in Lai.6-i,Ndo.4Sra;Cu04E3 
(x), ^^Cu wipeout onset Tnqr in Lai.6-a;Ndo.4Sra;Cu04 (a) 
and Lai.8-iEuo.2Sra:Cu04 (T), onset temperature T^hf 
for short . rMig e charge order according to neutron 
and X- rJBm in Lai.6-iNdo.4Sra;Cu04 (♦), the cop- 
per wipeout inflection point (■) in Lai.6-2,Ndo.4SriCu04, 
the temperature where = 0.05sec"^K" 

Lai.8-2,Euo.2Sr2,Cu04 (•), long range spin order Tspi, 
in Lai.6-a;Ndo.4Sra;Cu04 (o), onset of La wipeout for 
Lai.8|-:,Euo.2Sr2:Cu04 (A), onset of /iSR coherent precession 
r,iSjiE3 for Lai.6-xNdo.4Sr:cCu04 (o), the onset of "^La re- 
covery of signal in Lai.8-2,Euo.2Sr2,Cu04 (V), and the onset 
of ^^'®^Cu recovery of signal in Lai.8-2,Euo.2Sr2,Cu04 (□). 
The darker grey tones indicate increasingly slow fluctuation 
timescales. We also show the superconducting boundary as a 
white line and a dotted line that connects T^harge and the Cu 
wipeout inflection points as a guide for the eye. 



B. Results for x > | 



We have already discussed the sharp onset of the slow- 
ing of spin fluctuations that occurs at Tcharge for a; = 0.12 
which is caused by charge order. Looking at the higher 
hole concentrations, x = 0.16 and 0.20, we see similar be- 
havior in the 63 Cu NQR wipeout (Fig. [ish, " ^La NQR 
wipeout (Fig. and in (FigTy), although 

in each case the temperature scale is shifted increasingly 
downwards for the samples farther from the magic hole 
concentration a; ~ |. This naturally leads to a similar 
distribution of F that is shifted to lower temperature, 
as shown in Fig. [ic| c and [To|d. In the case of x = 0.16 
we find that all of the anchoring points of the distribu- 
tion are shifted to lower temperatures by ^^lOK while for 
X = 0.20 the shift is --20-30 K. The striking similarity of 
the shapes of F(r) for 0.12 < a; < 0.20 indicates that 



the fundamental physics of these systems is very simi- 
lar. In each of these cases we consider charge order as 
the cause of anomalous glassy slowing of spin fluctua- 
tions. With increasing hole concentration the tendency 
to charge order diminishes, causing the shift of F to lower 
temperatures for larger x. 

The same strong similarities for the samples with 
X > i can also be seen in the phase diagram of 
Fig. |l^ The lightest grey corresponds to the region 
where some segments of the Cu02 plane experience spin 
fluctuations slower than 63pu;zpeo«t^ 63 ^q^l sig- 
nal begins to wipe out. For Lai.6-xNdo.4Sr2;Cu04 with 
0.12 ^ X < 0.15, the boundary agrees well with the 
onset of short-range charge order as observed by scat- 
tering methods.LlO We note that our wipeout studies 
of Lai.6-£!;Ndo.4Sra;Cu04 were conductedtll without the 
knowledge of the more recent results from x-ray scatter- 
ing by NiemoUer et al. for x — 0.15,E3 and the good agree- 
ment between our Tnqh — 60±5 K and their Tcharge — 
55 K strongly supports our indentiflcation of T^qh as the 
onset of the glassy slowing of the spin fluctuations trig- 
gered by slowing charge dynamics for 0.12 < x < 0.15. 
Generally speaking, for x > ^, we find that all of the 
contours are nearly parallel and slope downward with 
increasing x. Thus below the onset temperature for 
charge order Tcharge, the slowing down is similarly rapid 
throughout this concentration range, but with an energy 
scale that decreases away from the magic doping le'we]. 
X ~ I where the stripe fluctuations are most robust .til 
In fact, if we fit the slowing of spin fiuctuations to the 
renormalized classical form exp(— 27r/9f-^-''/T) for different 
hole concentrations we find that decreases with 



increasing x above 



as shown in Fig. 19. 



C. Results for x < | 



Within the second region a; < ^, the slowing down 
is more gradual and the contours in Fig. |l^ are more 
spread out. At this point it is worth pointing out that 
although Tnqh increases below a; « |, ^^F{T) is more 
tailed as shown in Fig. |l^ for Lai.6-xNdo.4oSra;Cu04 with 
X = 0.07 and 0.09. It is important to note that there 
is an infiection point in the temperature dependence of 
the wipeout fraction, ^^F{T), at approximately 55 K and 
62 K, respectively, for these two concentrations. In a 
recent study on a Lai.5oNdo.4Sro.ioCuO^-&ingle crystal 
provided by Ichikawa etal, we confirmedt23 that this in- 
fiection point is nothing but the charge ordering teip=. 
perature Tcharge as determined by x-ray scattering.tJ 
The coincidence of Tcharge and the LTLO (Low Tem- 
perature Leas-|Orthorhombic) structural phase transition 
temperatureE3 may suggest that static one dimensional 
charge order cannot develop without the tetragonal lat- 
tice symmetry that is compatible with the symmetry of 
charge stripes. The longer spatial coherence of charge 
order below Tcharge accelerates wipeout at the inflection 
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point. Even though static charge order with substantial 
spatial coherence has not developed between Tjsiqr and 
Tcharge, transport measurements by Ichikawa etal. show 
that the onset of wipeout agrees well with the tempera- 
ture scale T„ for local charge order as deduced from a scaL 
ing analysis of the enhancement of in-plane resisitivity.LJ 
These results underscore the importance of the interplay 
between charge and spin degrees of freedom, even for 

X < I. 



ture. We include data for 27rp^^-^ for the range x < ^ 
in Fig. |l^, but we reiterate the differences between the 
physics above and below a; ~ ^ and use dotted lines 
for the underdoped samples. We note that the trend of 
It^pI^^ decreasing away from x ~ g is expected from the 
low temperature behavior seen in Fig. [l^, indicating that 
the relevant energy scale is indeed decreasing, even if an 
identical analysis is not entirely applicable. 
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FIG. 19. Doping dependence of 2-K(fJ^ as deduced from 
the data presented above in Fig. ^ by fitting to the form 
r oc exp( — 271/9^-'^ •'^/r). We note that F levels for low temper- 
ature and that region is neglected in the fit, resulting in large 
uncertainty in the value of li^p'^J^ . For a; < | tailed wipeout 
indicates that the physics may differ substantially from the 
region x > i, so we use dashed lines in that region. 

It would thus appear that for x < i, there is a pre- 
cursive mechanism that wipes out the ®'^Cu resonance 
signal below Tj^q^ before collective short range charge 
order sets in at Tcharge{< T^qr)- The nature of this 
precursive wipeout cannot be collective since there is no 
sharp onset, and th us c annot be explained in the,terms 
discussed in section IV, In earlier publications,E3E2l we 
have suggested that this precursive wipeout may be ac- 
counted for if we assume that individual holes begin 
to localize with decreasing temeprature, thereby creat- 
ing local moments. This is a situation similar to the 
one tkajb results in wipecoiijiiji conventional dilute Kondo 
alloysLj and spin glassesEJ'Lj On the other hand, in this 
scenario in which the local moments are created by lo- 
calization effects, the precursive wipeout may be further 
enhanced if one increases resistivity by introducing dis- 
order. However, our recent experimental studies on Zn 
doped La2-a;Sra;Cu04 show that the onset of wipeout 
Tnqr is a fairly well-defined temperature scale that does 
not vary even when resistivity is enhanced and-local mo- 
ments are introduced by Zn impurity doping.1123 Eurther- 
more, a recent inelastic neutron scattering studyll23 on a 
sample of Lai.93Sro.o7Cu04 showed a continual decrease 
of the fluctuation energy scale with decreasing tempera- 
ture, followed by the freezing of spin stripes at low tem- 
peratures. The presence of a striped ground state indi- 
cates that it is unlikely that the glassy spin fluctuations 
below T]\jQji have no spatially correlated stripe signa- 



VI. CONCLUSIONS 

In this paper, we provide a comprehensive NQR 
picture of the spatial modulation and fluctuations of 
spin and charge density waves in the stripe phase of 
Lai.88-y(Nd,Eu)j;Sro.i2Cu04 and Lai.875Bao.i25Cu04. 
From the analysis of the Zeeman perturbed ^^Cu and 
"'^^^La NQR spectra, we deduce the spatial distribution 
of the hyperfine field at the time scales of NQR spin 
echo measurements. We demonstrate that the maximum 
frozen Cu magnetic moment at 350 mK is relatively small, 
~ 0.1 — 0.2/iB, with a comparably large distribution. The 
fact that fiSK studies find a large and well defined effec- 
tive moment of '^O.SfiB at a timescale that is two or- 
ders of magnitude faster than NQR, provides evidence 
that motional avergaing is at least partially responsible 
for the supression of the effective moment as observed 
in our NQR measurements. Given that the mobility of 
holes in the Cu02 planea,at T ^ Tcharge is compara- 
ble to the metallic state,EJ perhaps it is reasonable that 
the magnitude of the frozen moment averaged over NQR 
time scales is fairly small and highly disordered, despite 
the fact that the spin-spin correlation length is known 
to reach ~ 200 A below 30 K at neutron scattering time 
scalesJj We also find that broadening of the quadrupole 
coupling is necessary to fit the low temperature Zeeman 
perturbed NQR linesape, indicating that there are spa- 
tial variations in the local charge state. 

From a simple physical argument based on comparison 
of and the renormalized classical scaling of the 

0(3) Non-Linear-cr model, we show that the roughly ex- 
ponential divergence of observed below Tcharge is 
consistent with the presence of low frequency spin fluctu- 
ations characterized by a reduced effective spin stiffness 
27rp'^^f Ri 200 K. We attribute the renormalization of the 
spin stiffness from 27rps ~ 10'^ K at higher temperatures 
to 2npl^f ~ 10^ K to slowing charge dynamics that sup- 
press the Cu-Cu exchange interaction. By analyzing the 
magnitude and distribution of the nuclear spin-lattice 
relaxation rate based the Non-Linear-cr model 
with 27Tp'^^f « 200 K, we deduce the spatial distribu- 
tion of F (spin fluctuation frequency scale) of Cu below 
Tcharge- Wc also analyze the ^^Cu and ^^^La NQR wipe- 
out fraction F{T), elastic neutron scattering, and fiSK 
data within the same framework, and obtain a nearly 
identical spectrum. Combined with the fact that only a 
small change of a factor of ^ 3 in the anomalous slowing 
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of r is sufficient to account for tfie observed wipeout of 
^^Cu NQR, tfiis provides a natural explanation of why 
the temperature dependence of F{T) tracks the charge 
order parameter measured by scattering techniques for 
0.12 < X < 0.16, which may be interpreted as the volume 
fraction of patches of the sample in which the spin fluc- 
tuations have started to slow exponentially. However, we 
note that our model is quite simple minded, relying en- 
tirely on magnetic enhancement of 1/Ti.t23 Fortunately, 
since the qualitative change in the temperature depen- 
dence of spin fluctuations from roughly T-linear to expo- 
nential is so drastic, the overall picture of the tempera- 
ture dependence of the spin fluctuation frequency scale 
F shown in Fig. ^ would not be substantially altered by 
a more elaborate analysis. 

By plotting contours of spin fluctuation time scales in 
the T — X phase diagram, we demonstrate that similar 
glassy slowing of stripes takes place away from the magic 
hole concentration x = However, we find a qualita- 
tive difference in the slowing between above and below 
X = i wheie the measured incommensurability begins 
to saturate JIj The origin of the difference is not well un- 
derstood at this point. For x < ^, stronger tendency 
toward some sort of localized stgytepseems to precede spa- 
tially periodic charge ordering j^^rn probably inducing 
local spin gla ss . li| ke behavior with fairly high tempera- 
ture scales J13'tll£3 However, we emphasize that this spa- 
tially disordered spin glass like entity eventually|-|freezes 
as stripes, as demonstrated by neutron scattering.la Hence 
Tnqr may be considered as the onset of the glassy freez- 
ing of stripes. A recent scaling analysis of resistivity 
data by Ichikawa et alt3 suggests that this glassy freez- 
ing is triggered by slowing charge dynamics, in support 
of our earlier proppsal for the mechanism of wipeout 
and glassy freezingO We also believe that the qualita- 
tive difference above and below a; « ^ bears significant 
implications on the mechanism of superconductivity be- 
cause the behavior of La2-a;Sr2;Cu04 is very similar to 
Lai.88-y(Nd,Eu)j,Sro.i2Cu04 at and below x — ^, while 
the glassy freezing of stripes is not seen by NQR in su- 
perconducting La2-2:Sr2:Cu04 for a; > |. 

Finally, we emphasize that the so-called stripe phase 
never looks striped at the slow time scales of NQR mea- 



surements even for 



at temperatures as low as 



350 mK. It is certaintly possible to argue that the seg- 
ments of the Cu02 plane in which the ^'^Cu NQR signal 
is wiped out are very similar to a conventional spin glass, 
and at slow time scales those segments are phase sepa- 
rated from other patches of the Cu02 plane where the 
Cu spins still behave as they did above Tch 



arge 



Whether 

such segments are fixed in space is not clear. The essence 
of the stripe physics is that these disordered patches of 
the Cu02 planes with slow spin and charge dynamics ap- 
pear striped if viewed by an experimental probe with a 
timescale that is faster than the remnant fluctuations. 
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